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I.  INTRODUCTION 


The  Interior  Ballistics  of  High  Velocity  Guns  (IBHVG2)  prograin  has  a  wide  and  loyal  following. 
U  has  developed  in  capabilities  and  sophistication  over  three  decades,  which  is  a  very  long  time  in 
computing.  IBHVG2  is  the  standard  against  which  many  other  interior  ballistics  modeling  projects  must 
first  demonstrate  their  accuracy.  In  view  of  this,  an  update  of  the  code  documentation  has  been  needed, 
since  the  "User’s  Guide"  (Fickie  and  Anderson  1987)  is  based  on  Version  4.00,  which  is  in  English  units, 
while  Version  5.C1,  which  is  presently  the  most  commonly  used  one,  is  in  metric  units. 

In  addition,  recent  propulsion  systems,  referred  to  as  electrothermal-chemical  (FTC)  guns  and 
generically  shown  in  Figure  1,  have  combined  electrical  energy  in  the  form  of  a  plasma  with  propellants. 
Endothermic,  exothermic,  solid,  liquid,  gel,  and  slurries  have  been  considered  for  the  oropellanL  Although 
various  designs  have  been  tested,  it  is  useful  to  develop  a  generic,  zero-dimensional  simulation  capability 
to  explore  both  the  physics  ~>f  the  plasma-propellant  interaction  and  baseline  performance  predictions. 


Figure  1.  Generic  elcctrothennal-chemical  (ETC)  gun. 

IBHVG2  is  a  lumped  parameter  (zero-dimensional)  interior  ballistic  code  which  has  been  distributed 
by  the  Weapons  Technology  Directorate,  Army  Research  Laboratory  (WTD/ARL),  formerly  the  Ballistic 
Research  Laboratory  (BRL),  to  researchers  involved  with  gun  propulsion  technology.  The  code  is  capable 
of  modeling  the  ballistic  cycle  of  multiple  solid  propellants  from  propellant  ignition  to  the  blowdown  of 
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combustion  gases  after  the  projectile  exits  from  the  muzzle.  Various  automated  search  techniques  are 
available  for  finding  optimized  gun  parameters,  desired  propellant  characteristics,  or  matching  experimental 
data.  The  code  also  has  an  expansive  output  capability  which  allows  the  user  to  analyze  variables  of 
interest  over  the  time  of  the  simulation.  All  the  available  options  are  activated  from  an  input  file,  making 
IBHVG2  an  ideal  tool  for  the  scientist  who  does  not  wish  to  program. 

This  tutorial  provides  current  information  on  how  to  use  the  metric  version  of  IBHVG2.  The  first  part 
of  the  tutorial  describes  the  key  differences  between  Version  5.04  and  Version  4.00.  It  also  includes  a 
description  of  extensions  to  the  code  for  electrothermal  (ET)  and  electrothermal-chemical  (ETC) 
calculations.  The  second  part  of  the  tutorial  provides  an  updated  summary  of  the  input  decks  and  presents 
sample  code  calculations  using  the  ETC  enhancements.  The  major  enhancements  are  (1)  revision  of  the 
energy  equation  to  account  for  electrical  energy;  (2)  provision  for  description  of  input  power,  and  (3) 
provision  for  a  generalized  gas  generation  rate.  All  search  techniques  in  the  original  code  have  been 
maintained. 

2.  THE  INPUT  RLE 

The  input  format  for  IBHVG2  consists  of  an  input  file  composed  of  various  descriptive  inpit  decks. 
All  the  input  deck  titles  and  variables  are  written  in  uppercase  letters.  Each  input  deck  is  qualified  by  a 
dollar  sign  in  column  1 ,  followed  immediately  by  the  title  of  the  deck.  A  dollar  sign  in  any  other  column 
indicates  to  the  code  that  what  follows  on  that  line  is  a  comment.  No  tabs  are  allowed  anywhere  in  the 
input  file. 

Within  an  input  deck,  variables  are  defined  by  an  equals  sign  followed  by  the  value,  with  spaces  about 
the  equals  sign  squeezed  out.  Some  variable  names  have  synonyms  which  can  be  found  in  the  "User’s 
Guide."  Arrayed  variables  can  be  allocated  values  either  by  individually  defining  them. 

GL0C(1)=10 

GLOC(2)=20 

GLOC(3)=30 


or  by  separating  the  values  by  commas, 


GLOC  =  10,  20.  30 
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noting  that  a  space  after  a  value  (instead  of  a  comma)  indicates  the  end  of  the  variable’s  definition.  When 
multiple  entries  are  contained  within  an  input  deck,  the  code  uses  the  last  definition  for  the  simulation. 

Character  strings  must  be  included  in  quotes  or  apostrophes  when  defined,  for  example, 

SHOW  =  *GAGE(2)' 

The  most  recent  version  of  IBHVG2  distributed  by  the  ARL  is  version  5.04,  which  uses  metric  units. 
The  standard  measure  of  length  is  in  meters,  weight  is  in  kilograms,  and  pressure  is  in  MegaPascals.  All 
input  and  output  parameters  are  described  in  these  units.  The  only  exception  to  this  is  when  a  user 
references  a  pressure  value  from  a  STDIS  or  $PD1S  deck.  In  these  cases,  the  referenced  pressure  is  in 
Pascals. 

A  summary  of  the  input  decks  available  in  an  input  file  follows. 

REQUIRED  INPUT  DECKS 

$GUN  -  Defines  gun  physical  dimensions  and  pressure  probe  locations. 

SINFO  -  Defines  error  tolerance  of  run,  pressure  gradient  selection,  times tep  of  simulation,  and  output 
print  options. 

SPROJ  -  Defines  projectile  parameters  (mass). 

SPRIM  --  Defines  primer,  which  is  assumed  all  burnt  at  lime  0. 

SPROP  --  Defines  main  propellant  charge(s)  (grain  dimensions,  physical  properties,  bum  rate 
and  ijr ution  behavior). 

SEND  -  Defines  end  of  input  file. 

OPTIONAL  INPUT  DECKS 

SRESI  --  Defines  gun  tube  resistance  to  projectile  motion  (shot  start  resistance),  and  air  resistance  to 
projectile  in  gun  tube. 

SHEAT  --  Defines  heat  loss  from  combustion  gases  to  a  thin  shell  of  the  gun  tube  (included  by  default). 
$ETC  -  Defines  electrical  power  input,  and,  if  desired,  lime-dependent  gas  generation  rate. 

SRECO  -  Defines  weight  of  gun  which  is  used  to  calculate  an  energy  loss  due  to  recoil. 

SCOMM  -  Comment  deck. 

SSAVE  -  Used  after  the  SEND  card,  saves  previous  input,  and  allows  the  user  •."*  change  only  desired 
parameters  for  another  run. 
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PRINTOUT  OPTIONS 


STDIS  --  Specifies  user-selected  trajectory  variables  for  printing. 

SPDIS  -  Specifies  user-selected  variables  to  be  printed  after  each  parametric  run. 

SEARCH  OPTIONS 

SFTND  -  Searches  on  any  parameter  to  give  a  desired  output  variable. 

$PARA  -  Performs  a  systematic  variation  of  any  parameter. 

SPMAX--  Searches  on  a  desired  propellant  characteristic  to  give  a  selected  maximum  breech  pressure. 

The  following  are  not  input  decks,  but  are  lists  of  variables  to  choose  from  for  output  printing  and 
for  the  SFIND  search  option. 

STRAJ  -  List  of  trajectory  variables  (see  Appendix  A).  Additional  variables  are  available  to  include 
ETC  parameters,  see  $ETC  deck  description. 

SOUT  -  List  of  output  variables  (see  Appendix  A). 

3.  DESCRIPTION  AND  DISCUSSION  OF  THE  INPUT  DECKS 

3.1  SGUN  Deck.  The  SCUN  deck  describes  the  gun  tube  and  chamber  physical  dimensions  to  be 
used  by  the  code.  The  chamber  volume  (CHAM),  tube  diameter,  and  length  must  be  specified.  If  the  gun 
tube  is  smooth-bore,  the  diameter  is  defined  by  the  LAND  parameter.  (See  Figure  1.)  The  TRAV 
parameter  defines  the  length  of  the  gun  barrel  from  the  rear  sealing  point  of  the  projectile  to  the  muzzle. 

If  a  gun  tube  with  rifling  is  to  be  modeled,  the  diameter  of  the  groove  (GRVE)  must  be  input  with 
a  value  of  twist  (TWST)  in  calibers  per  turn.  A  value  for  the  groove  over  land  surface  area  (GIL)  must 
also  be  included.  (See  Figure  2.)  The  code  uses  these  parameters  to  calculate  a  rotational  kinetic  energy 
of  the  projectile. 

Pressure  probe  locations  in  the  gun  tube  and  chamber  are  specified  in  the  SGUN  deck.  NGAG  defines 
:he  number  of  pressure  probes  (up  to  30),  and  GLOC  is  an  array  of  size  NGAG  which  defines  the 
locations  of  each  respective  probe.  Distances  are  measured  from  the  projectile  base.  into  the  chamber 
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GUN  TUBE  STEEL 


LAND 


Figure  2.  Additional  parameters  required  for  »  gun  with  a  rifled  bore- 

or  "+"  down  the  barrel.  The  effective  length  of  the  chamber  (CLEN).  should  be  input  to  the  code  when 
using  in-chamber  pressure  probes.  For  example.  Figure  3  contains  the  simulated  pressure  probe  locations 
used  for  30-mm  experiments  at  GT-Devices  (McElroy,  Grieg,  and  Juhasz  1991). 

The  maximum  pressure  calculated  by  the  code  in  the  gun  at  any  given  time  is  always  located  at  the 
breech  face.  The  pressure  at  the  in-duunber  pressure  probe  locations  are  scaled  down  using  the  selected 
pressure  gradient  model  (chosen  in  the  SINFO  deck).  The  code  calculates  the  pressure  scaling  for  the  in¬ 
chamber  pressure  probes  using  CLEN,  which  by  default  is  the  chamber  volume  divided  by  the  bore  area. 

If  chambrage  is  present  in  the  gun  chamber,  1BHVG2  has  been  modified  to  calculate  the  effect  of  the 
change  of  diameter  in  the  shoulder  region  on  the  developed  pressure  (Robbins,  Gough,  and  Anderson 
1989).  The  chambrage  pressure  gradient  model  is  selected  in  the  SINFO  deck,  but  the  physical 
dimensions  are  entered  in  the  SGUN  deck.  The  location  of  the  diameter  change,  DIST,  is  referenced  from 
the  breech  face  The  diameter  at  the  location  corresponding  to  DIST  is  input  as  DIAM.  The  number  of 
location-diameter  pairs  is  input  as  CPTS,  with  a  maximum  of  five  allowed.  (See  Figure  4.) 
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Number  of  Pressure  Proto**  -  NOAQ  ■  6 

Position* -OLOC.-O.1ie,  -0JDS5,  OJM,  0.148,  0-334,  1.301 
PCI  PC2  PS1  PB2  PB3  PB4 


3.2  S1NFO  Deck.  The  SINFO  deck  defines  the  run  options  for  a  simulation.  The  title  of  the  run  (up 
to  48  characters)  is  defined  in  this  input  deck  by  the  variable,  RUN.  and  is  printed  at  the  top  of  each  page 
of  the  output. 

The  maximum  integration  timestep  is  set  by  DELT,  and  the  printout  timestep  is  set  by  DELP. 
Generally,  if  comparison  to  experimental  data  is  desired,  it  is  useful  to  set  the  printout  timestep  to  match 
the  digitizer  sampling  rate. 

If  printout  is  desired  in  steps  of  projectile  travel,  set  the  parameter  DLPU  equal  to  2  and  note  that 
DELP  is  now  the  printout  step  in  units  of  projectile  travel. 

EPS  is  the  maximum  allowed  error  for  the  integrator  timestep  adjustment  in  Runga-Kutta  integration 
subroutine.  This  parameter  ultimately  controls  the  run  time  for  convergence  of  the  program.  If  the  code 
terminates  during  a  simulation  or  continues  to  run  for  an  excessively  long  time,  the  EPS  parameter  is  the 
most  suspect  Generally,  if  the  EPS  parameter  is  increased,  the  code  converges  to  a  solution  quicker. 
Conversely,  if  greater  accuracy  is  desired,  this  parameter  should  be  decreased. 

The  models  used  for  determining  the  pressure  drop  down  the  gun  tube  from  the  breech  face  can  be 
chosen  from  the  Lagrange  (GRAD=I).  Pidduck-Kent  (GRAD=2),  or  chambrage  (GRAD=3)  pressure 
gradients.  The  Lagrange  pressure  gradient  is  appropriate  for  propellant  charge  over  projectile  mass  (C/M) 
ratios  of  less  than  2.  The  Pidduck-Kent  pressure  gradient  is  often  utilized  for  C/M  ratios  of  greater  than 
one  and  is  the  best  choice  for  large  C/M  ratios,  if  chambrage  is  present,  a  pressure  gradient  model 
incorporating  the  changes  in  chamber  diameter  is  available.  The  physical  dimensions  of  the  diameter 
changes  in  the  chamber  is  input  in  (he  5GUN  deck. 

Output  file  print  options  are  specified  by  the  POPT  variable  array,  which  contains  five  elements.  The 
default  values  for  each  element  is  1,  except  POPT(4)  which  has  a  default  of  0. 

POPT(I)  set  equal  to  1  echoes  the  input  file  specifications  at  the  beginning  of  the  output  file  (setting 
to  0  turns  off  this  option).  This  can  be  of  help  in  locating  errors  in  the  input  format,  since  the  code  will 
locate  by  card  number  and  column  an  input  it  docs  not  understand. 
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The  second  element  of  the  array.  POPT(2),  represents  the  type  of  "trajectory''  variables  vs.  time  a  user 
might  want  to  specify  for  a  simulation.  POPT(2)  set  equal  to  0  instructs  the  code  to  skip  the  trajectory 
printout.  When  POPT(2)  is  set  equal  to  1,  trajectory  print  includes  time,  projectile  travel,  projectile 
velocity,  projectile  acceleration,  breech  pressure,  mean  pressure,  base  pressure,  mean  temperature,  and  the 
fraction  of  each  propellant  deck  burned.  The  user  can  select  to  print  out  desired  trajectory  variables  from 
the  STRAJ  list  by  setting  POPT(2)  equal  to  2.  Each  user-defined  output  trajectory  variable  must  then  be 
specified  in  a  STDIS  deck,  which  allows  a  multiplication  or  division  of  the  variable  to  accommodate  a 
change  in  units.  Up  to  11  variables  are  allowed  to  be  specified  via  STDIS  specifications. 

The  third  variable.  POPT(3),  when  set  equal  to  1,  activates  a  summary  of  the  simulation  which 
provides  the  default  trajectory  values  at  both  the  maximum  breech  pressure  and  at  muzzle  exit.  This 
summary  is  very  convenient  since  it  provides  ballistic  performance  in  compact  form.  Setting  equal  to  0 
turns  off  this  option. 

POPT(4)  provides  the  user  a  choice  of  including  a  pressure  blowdown  phase  after  the  projectile  leaves 
the  muzzle.  Set  equal  to  0,  the  code  ignores  this  option.  When  POPT(4)  equals  1.  the  code  calculates 
blowdown  until  the  rarefaction  wave  reaches  the  breech  face.  When  POPT(4)  equals  2  or  3.  an  additional 
parameter.  BLPR.  must  be  included.  This  parameter  limits  the  blowdown  calculation  to  a  desired  breech 
pressure.  The  results  of  the  blowdown  calculations  are  printed  at  the  end  of  the  output  file. 

POPT(5)  is  a  print  option  similar  to  POPT(2),  except  that  it  provides  output  selections  for  parametric 
variation  runs.  When  POPT(5)  equals  0.  the  previously  defined  trajectory  printout  options,  via  POPT(2) 
specifications,  are  used  for  each  result  of  the  parametric  calculations.  When  POPT(5)  is  set  equal  to  1. 
the  code  prints  a  full  trajectory  output  for  the  first  run  and  a  single  line  summary  for  each  additional  run. 
POPT(5)  set  equal  to  2  activates  a  user-defined  parametric  summary  printout  for  each  run  via  SPDIS 
specifications.  Each  SPDIS  deck  defines  one  variable  to  be  printed.  SPDIS  decks  are  similar  to  STDIS 
decks,  although  the  variable  selection  is  made  from  input  variables.  STRAJ  list,  or  in  the  case  of  the  single 
line  summary,  the  SOUT  list.  The  output  selections,  for  example,  can  be  made  to  contain  the  changing 
input  parameters  as  well  as  desired  output  values.  A  summary  of  the  print  options  is  contained  in 
Figure  5. 
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POPT(l) 

=  0 

No  input  echo 

=  1 

Echoes  input 

POPT(2) 

=  0 

No  trajectory  printout 

=  1 

Default  trajectory  printout 

=  2 

User-defined  trajectory  printout  via  STDIS  specifications 

POPT(3) 

a  0 

No  interior  ballistics  summary 

a  1 

bterior  ballistics  summary 

POPT(4) 

=  0 

No  blowdown  calculations 

a  t 

Blowdown  calculations  until  rarefaction  wave  reaches  breech  face 

a  2 

Blowdown  calculations  until  reaching  a  desired  breech  pressure,  BLPR 
(reduced  printing 

=  3 

Blowdown  calculations  until  reaching  a  desired  breech  pressure,  BLPR 
(expanded  printing) 

POPT(5) 

=  0 

Printont  for  each  parametric  calculations  in  format  consistent  with 
POPT(2)  specification 

a  1 

Printout  for  parametric  calculations  with  a  complete  trajectory  printout 
for  the  first  ran  and  single  line  summary  for  each  additional  parametric 

nm 

a  2 

run 

User-defined  parametric  summary  printout  via  SPDIS  specifications 

FiguA.  5.  Summary  of  POPT  options. 


The  in^nr.*  >.  n  folio**  inc  fte  input  echo  and  before  the  trajectory  printout  is  the  codes  interpretation 
of  the  variabi-s  aefimv;  ;n  the  iuput  file.  These  specifications  represem  the  values  the  code  will  use  for 
the  simulation.  The  code  is  very  capable  (through  default  values,  user  input  errors,  etc.)  of  misinterpreting 
a  user  s  intentions  without  generating  a  parsing  error. 


A  table  of  the  energy  distribution  of  the  simulation  is  contained  at  the  end  of  the  output  file  after  the 
internal  ballistics  sumnwy.  Various  parameter*  are  calculated  after  the  run  and  are  printed  at  the  end  of 
the  input  file:  loading  density.  C/M,  piezometric  efficiency  (mean  pressure/max.  breech  pressure), 
electrical  enhancement  factonproj.  KE/electrical  energy),  v*d  the  expansion  ratio  (chamber  volume  +  bore 
volumc/chamber  volume). 


To  obtain  the  maximum  efficiency  of  a  specific  gun  fixture,  a  constant  pressure  option  is  also  available 
in  this  input  deck.  The  code  holds  a  constant  pressure  during  the  ballistic  cycle  by  varying  the  burning 
rate  (CONP»l)  or  the  surface  area  (CONP-7.)  of  a  single  propellant  The  desired  constant  pressure  is 
defined  by  PRES,  with  an  error  tolerance  on  the  pressure  defined  as  TOL.  Only  one  propellant  deck 
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containing  a  homogeneous  (non-layered)  grain  can  be  included  in  this  simulation.  The  propellant  deck 
must  begin  burning  at  the  start  of  the  simulation  with  no  ignition  delay. 

3.3  SPROJ  Deck.  This  deck  defines  projectile  parameters  for  the  simulation.  The  only  required 
parameter  is  the  projectile  weight,  PRWT.  Parameters  for  complex  projectiles  are  also  available,  but  will 
not  be  discussed. 

By  specifying  the  C/M  parameter,  an  option  of  calculating  either  the  projectile  weight  given  a  total 
charge  weight  (CO)PT=l).  or  calculating  a  total  charge  weight  given  a  projectile  weight  (COPT=2)  is 
available.  If  COPT=2  is  used,  the  calculated  total  charge  weight  is  distributed  according  to  the  relative 
percentages  of  the  initial  individual  charge  weights. 

3.4  iPRIM  Deck.  The  primer  deck  is  required  by  the  code  for  each  simulation.  Required  parameters 
are  the  specific  heu  "ado  (GAMA),  impetus  (FORC),  covolume  (COV),  flame  temperature  (TEMP),  and 
charge  weight  (CHWT).  Enough  energetic  material  should  be  included  to  raise  the  chamber  pressure  to 
at  least  0.34  MPa  (50  psi).  The  primer  is  assumed  to  be  completely  burned  at  time  0. 

3.5  SPROP  Deck.  The  propellant  deck  defines  the  main  propellant  charge  characteristics.  Up  to  five 
different  decks  are  allowed  by  the  code.  The  grain  density  (RHO),  GAMA.  FORC.  COV,  TEMP,  CHWT. 
grain  geometric  shape,  and  burning  rate  must  be  specified  for  each  charge. 

Various  form  functions  for  defining  propellant  grain  geometries  (GRAN)  can  be  chosen,  including  a 
completely  general  form  function  which  is  defined  by  surface  areas  at  a  depth  burned.  The  simplest  form 
function  is  the  spherical  grain  (GRAN«"SPHR"),  which  only  requires  a  value  of  D1AM  for  physical 
description.  The  form  function  for  ball  powder  grains,  such  as  WC891,  is  GRAN»"CAKE",  which 
requires  a  DIAM  and  a  thickness  (THCK)  for  physical  description.  The  perforated  grains  used  in  large- 
caliber  guns  for  increased  surface  area  require  the  most  detail  (inner  and  outer  web,  perforation  diameter, 
outer  dimensions,  etc.)  for  complete  physical  description.  Fortunately,  examples  and  illustrations  for  the 
more  complex  grain  geometries  are  contained  in  the  "User’s  Guide.”  The  general  form  function  is 
specified  by  the  number  of  depth  burned-surface  area  pairs  (NSUR).  with  the  array  of  depths  defined  as 
DEPB  and  the  array  of  surface  areas  at  the  corresponding  depth  burned  defined  as  SURF.  Ratios  of  a 
grain  dimension  to  one  or  the  other  given  dimensions  are  available  which  allow  optimization  via  the 
search  options  to  find  desired  propellant  grain  size  characteristics. 
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The  code  in  its  present  form  can  describe  a  variation  of  propellant  properties  (thermochemistry  and 
bum  rate)  within  a  single  grain  with  up  to  four  points  at  desired  depths  into  the  grain.  Descriptive 
information  is  input  for  each  layer  boundary  into  arrayed  variables  of  length  four,  which  have  different 
names  than  those  previously  discussed  for  homogeneous  grains.  Density  at  a  layer  boundary  is  now  input 
as  RHOL,  specific  heat  ratio  as  GAML,  impetus  as  FRO,,  covolume  as  COVL,  and  (lame  temperature 
as  TMPL.  An  example  for  a  four-layer  density  input  is  as  follows: 

RHOL  =  1000.  1200.  1300.  1350 

density:  layer  1.  layer  2.  layer  3.  core. 

The  last  letter  in  all  of  these  layered  variables  corresponds  to  the  surface  in  question:  L  for  lateral, 
E  for  end.  and  P  for  perforation.  In  th;  case  of  the  CAKE  form  function,  i.e.,  ball  powder,  only  lateral 
surface  definitions  are  required.  This  is  also  the  case  for  spherical  grains. 

Transition  depths  from  layer  to  layer  are  required  and  are  input  into  the  DEPL  array  of  length  3. 
These  transition  depths  from  one  layer  boundary  to  the  next  are  flagged  during  the  simulation  in  the 
trajectory  printout  as  the  bum  progresses  into  the  next  layer.  An  example  for  a  spherical  grain  is 
presented  in  Figure  6. 


Spherical  grain 

DEPI.  =  0.0001,  0.0002,  0.0003 
DIAM  =  0.002 

Figure  6.  Details  of  the  transition  depths  for  a  layered  propellant  grain. 
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For  propellants  which  can  be  defined  with  less  than  four  layers,  the  user  is  still  required  to  all  the  four 
clement  arrays  (for  the  propellant  properties),  although  it  is  sufficient  to  just  duplicate  the  values  for  the 
remaining  unused  layers.  If  only  one  layer  is  present,  it  is  taken  to  be  the  core  layer  and  only  one  set  of 
values  js  needed.  The  transition  depth  array  must  contain  values  of  all  three  elements  which  are  distinct 
and  increasing. 

Note  that  the  code  linearly  interpolates  the  given  values  from  one  layer  boundary  to  the  next,  including 
bum  rate,  with  the  center  of  the  grain  Gayer  4)  having  constant  properties  throughout  An  example  of  the 
code  interpretation  for  the  variation  in  density  for  a  layered  grain  is  shown  in  Figure  7. 

Propellant  bum  rates  can  be  specified  either  by  direct  input  of  bum  rate  versus  ptessur  or  by  defining 
coefficients  and  exponents  for  the  classical  exponential  bum  rate  function  r  =  pp01. 


RHOL  s  1000, 1200, 1300, 1400 


Radius  of  grain 

Figure  7.  Example  of  the  interpolation  scheme  for  a  laveted  grain. 
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To  input  the  bum  rate  as  a  function  of  pressure,  set  NTBL  to  the  number  of  bum  rate-pressure  pairs. 
PR1L  is  the  pressure  array,  and  BRIL  is  the  bum  rate  at  that  respective  pressure. 

If  the  grain  is  homogeneous  and  the  bum  rate  is  desired  to  be  specified  in  the  exponential  form  by 
a  single  coefficient-exponent  pair,  the  coefficient  variable  is  BETA,  and  exponent  variable  is  ALPH 
(NTBL  can  be  ignored  or  set  equal  to  0). 

Otherwise,  to  input  the  bum  rate  in  the  exponential  form,  set  NTBL  to  the  number  of 
coefficient/exponent-pressure  pairs  but  make  it  negative.  The  minus  sign  signifies  to  the  code  that  the 
desired  entry  format  is  the  exponential  bum  rate  function.  PR1L  is  the  pressure  array,  CFIL  is  the 
coefficient  array,  and  EX1L  is  the  exponent  array  at  the  corresponding  pressure.  If  the  pressure  array  is 
deleted,  the  values  of  the  coefficient  and  exponent  are  constant  for  all  pressures  (NTBL=~1). 

The  numbers  in  these  variable  names  correspond  to  the  referenced  layer.  If  a  homogeneous  grain  is 
modeled,  the  code  requires  all  the  necessary  input  in  the  array  at  the  location  corresponding  to  the  fourth 
layer  boundary  (core).  Therefore,  the  number  in  these  variables  should  be  set  equal  to  4,  for  example, 
PR4L,  BR4L,  or  PR4L;  CF4L.  and  EX4L.  The  last  letter  corresponds  to  the  surface,  as  previously 
discussed. 

The  propellant  decks  by  default  will  begin  burning  at  time  0,  although  an  ignition  delay  can  be 
incorporated  into  a  simulation  through  a  variety  of  means.  The  variable  IGNC  in  the  propellant  deck 
controls  the  ignition  specification.  A  SPROPdeck  may  be  ignited  at  a  specified  time  (IGNC=I).  projectile 
position  (IGN02),  mean  pressure  (!GNC=3).  or  at  a  mass  fraction  burned  of  the  previous  propellant  deck 
(IGNC=4).  The  threshold  values  for  these  specifications  is  input  via  the  THRC  variable. 

An  additional  ignition  variation  miy  be  activated,  1GNS,  which  ignites  specific  surfaces  of  the 
individual  propellant  grains  at  thresholds  specified  by  THRS.  IGNS  and  THRS  are  arrays  of  length  3, 
with  each  element  corresponding  to  the  P,  E.  and  L  grain  surfaces,  respectively.  The  same  ignition  codes 
arc  used  as  with  IGNC  variable. 

The  propellant  deck  contains  many  options  and  is  therefore  the  most  complicated  input  deck.  Errors 
are  easily  made,  and  the  user  should  carefully  inspect  the  code  interpretation  of  the  input  values  in  the 
output  file,  especially  for  a  newly  created  propellant  deck.  The  user  will  find  propellant  decks  take  the 
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most  time  to  create  and  debug.  Therefore,  using  previously  created  ones,  if  available,  is  desirable. 
Sample  test  cases  containing  various  propellant  decks  can  be  found  at  the  end  of  the  "User's  Guide." 

3.6  SEND  Canj.  The  end  of  an  input  file  is  signified  to  IBHVG2  by  an  SEND  card.  No  additional 
lines  are  allowed  after  an  SEND  card,  except  a  SSAVE  deck. 

3.7  SRESI  Deck.  Gun  tube  resistance  to  the  motion  of  the  projectile  is  input  via  the  resistance  deck. 
The  resistance  profile  can  be  very  important  in  an  interior  ballistics  simulation  because  it  controls  how 
fast  the  combustion  volume  is  increasing.  This  affects  the  mean  pressure,  upon  which  the  propellant  bum 
rate  is  highly  dependent. 

NPTS  is  the  number  of  distance-resistance  pressure  pairs,  with  a  maximum  of  20  points  allowed. 
TRAV  (array  of  length  NPTS)  is  the  distance  from  shot  start  at  which  a  given  resistance  to  motion  occurs, 
and  PRES  (array  of  length  NPTS)  is  the  magnitude  of  the  resistance  pressure  at  the  respective  distance. 
Note  that  the  code  linearly  interpolates  between  the  resistance  pressure  points. 

A  multiplication  factor,  FACT,  is  available  to  quickly  scale  the  resistance  pressure  profile  up  or  down 
to  match  experimental  data.  RFPT  is  a  parameter  which  defines  the  first  point  in  the  resistance  array  to 
begin  multiplying  by  FACT. 

Air  resistance  to  the  motion  of  the  projectile  in  the  gun  tube  is  calculated  by  die  code  when  the 
parameter  AIR  is  set  equal  to  1 . 

3-8  SHEAT  Deck.  Heat  loss  from  the  combustion  gases  to  the  gun  barrel  is  defined  in  the  SHEA7 
input  deck.  The  parameter  HL  is  the  toggle  switch  for  including  heat  loss  (HL=1).  To  turn  off  this 
option,  set  HL=0.  Heat  is  transferred  from  the  combustion  gases  to  the  barrel  through  a  convective  heat 
transfer  coeffi:icni,  HO  (default  is  1 1.35  W/m2-  K).  Over  the  lime  scJc  of  >h  '’muiations  (milliseconds), 
the  gun  band  heat  sink  is  modeled  as  a  thin  shell  (default.  TSHL=0. 1  .  a  given  specific  heat. 

The  shell  is  initially  at  room  temperature  (TWAL=293  K)  by  default  For  guns  in  different  Uv*rre?l 
environments,  the  initial  temperature  can  be  defined  with  the  TWAI  parameter.  Default  vJues,  which 
cover  a  wide  range  of  steels,  are  given  for  the  gun  tube  density  (RSHL=7361.1  kg/mJ)  and  si'ecific  heat 
(CSHL=460.3  J/K). 
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Ur.'css  changes  arc  desired  for  die  input  parameters  given  above,  the  inclusion  of  the  JHEAT  input 
deck  in  an  input  file  is  not  required;  by  default  ail  the  calculations  are  performed. 

The  code  is  capable  of  modeling  only  one  heat  sink,  therefore  Lexan  chamber  liners  and  capillary 
liners  present  in  ET  guns  cannot  be  modeled  in  addition  to  the  steel  gun  tube. 

3.9  SETC  Deck.  In  the  solid  propellant  electrothermal -chemical  (SPETC)  gun,  electrical  energy  in 
the  form  of  a  plasma  is  utilized  in  the  combustion  chamber  to  perform  various  roles.  A  discussion  of  ETC 
guns  can  be  found  in  the  literature  (Juhasz  et  al). 


The  parameters  incorporated  into  this  deck  arc  electrical  power  addition  and  time-dependent  gas 
generation  rate  of  the  propellant.  Each  of  these  parameters  is  input  as  a  piecewise  linear  function  of  time. 
This  feature  allows  easy  input,  for  example,  from  an  experimental  gun  electrical  power-time  curve  by  just 
entering  points  from  the  experimental  power  curve. 

The  addition  of  electrical  power  input  directly  into  IBHVG2  updates  the  previous  method  which 
required  the  user  first  to  input  the  electrical  energy  in  the  form  of  electrons  into  the  BLAKE  (Freedman 
1982;  Bunte  and  Oberle  1989)  thermochemistry  code.  The  output  parameters  from  BLAKE  on  the 
electrically  augmented  propellant  were  then  input  into  IBHVG2.  This  laborious  process  is  no  longer 
needed,  assuming  the  electrical  energy  does  not  alter  propellent  thermochemistry  or  propellant  bum  rate, 
and  that  propellant  gas  proper  Jes  determined  at  a  fixed  loading  density  are  applicable.  However,  it  is 
noted  that  the  BLAKE  calculations  show  that  this  assumption  may  not  be  strictly  valid  (Wren  and  Oberle 
1992).  With  the  present  code,  ETC  gun  electrical  energy  augmentation  is  modeled  directly  in  IBHVG2 
with  unmodified  propellant  thermochemistry. 


The  electrical  power  input.  P,  is  provided  in  the  SETC  deck  a 3  a  piecewise  linear  function  of  time. 
This  function  is  integrated  in  closed  form,  and  the  resulting  energy  is  added  to  the  left  side  of  the  energy 
balance  equation,  Eq.  (A.01),  which  becomes 
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The  complexity  of  specifying  bum  rates  as  a  function  of  pressure,  or  providing  coefficients  and 
exponents  for  the  bum  rate  can  be  overridden  to  allow  input  of  a  gas  generation  rate  as  a  function  of  time. 
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An  exotic  propellant  can  be  modeled  without  attempting  to  determine  both  the  bum  rate  and  the  surface 
areas  by  providing  an  estimate  of  the  gas  generation  rate  versus  time.  The  system  designer  can  then 
perform  simulations  to  determine  a  gas  generation  rate  which  meets  system  requirements,  and  then 
subsequently  design  the  propellant  surface  areas  and  bum  rates  which  can  produce  the  optimal  gas 
generation  rate.  The  gas  generation  rate  is  also  specified  as  a  piecewise  linear  function  of  time.  This 
function  replaces  rh^  in  Eq.  (A.  15). 

The  input  parameters  related  to  ETC  gun  modeling  are: 

ELECTRICAL  POWER 

NPWR  -  integer  number  of  power-time  pairs  (min  =  2,  max  =  20) 

TPWR  -  time  (s) 

PWR  -  power  in  Watts  at  time  TPWR 

GAS  GENERATION  RATE 

NBRN  -  integer  number  of  gas  generation  rate-time  pairs  (min  =  2.  max  =  20) 

TBRN  -  time  (s) 

BRN  •  grs  generation  rate  at  lime  TBRN  (kg/s) 

The  following  related  parameters  have  been  added  to  the  STRAJ  deck: 

EPWR  -  Electrical  power  (Watts) 

EENE  -  Electrical  energy  integrated  from  t=0  (Joules) 

XBRN  -  Gas  generation  rate  (kg/s) 

Any  or  all  of  these  can  be  specified  in  a  STDIS  deck,  and  the  corresponding  variable  will  be  printed 
out  as  a  function  of  time  during  the  run. 

3.10  SRECO  Deck.  Recoil  losses  for  a  gun  system  can  be  calculated  in  the  code  through  the  use  of 
the  recoil  input  deck.  The  only  required  parameter  is  the  weight  of  the  recoiling  system 
(RCWD.Gcnerally.  the  energy  expended  in  recoil  of  a  gun  is  of  the  order  of  a  tenth  of  a  percent  of  the 
total  energy,  therefore  the  inclusion  of  recoil  in  a  simulation  has  a  negligible  effect. 


16 


3.11  SCOMM  Deck.  This  input  deck  is  used  for  making  comments  which  identify  a  simulation. 
Everything  that  follows  until  the  next  input  deck  is  considered  a  comment  Multiple  comment  decks  are 
allowed. 

3.12  SSAVE  Deck.  The  SSAVE  input  deck,  after  a  SEND  card,  allows  the  user  to  save  all  the 
previous  infonnation  given  in  the  input  decks  and  change  only  desired  variables  to  make  a  variation  on 
the  previous  simulation  all  in  the  same  tun.  The  input  deck(s)  which  contains  the  variable(s)  to  be 
changed  must  referenced,  followed  by  the  variable  redefinition.  Another  SEND  card  must  follow  the 
changes  to  indicate  to  the  code  that  the  end  of  the  input  has  been  reached. 

3.13  STD1S  Deck.  Each  STDIS  deck  defines  a  single  trajectory  variable,  from  the  STRAJ  list,  to  be 
printed  in  the  output  as  a  function  of  time.  A  maximum  of  1 1  STDIS  decks  are  allowed  in  a  simulation. 
The  parameter  SHOW  identifies  the  variable  to  be  printed.  A  multiplication  (MULT)  or  division  (DIV) 
of  the  selected  output  parameter  is  available  to  allow  a  change  in  units.  At  the  top  of  each  printed  page 
of  the  trajectory  output,  a  remark  (REM1)  is  printed  which  helps  identify  each  column.  An  example  is 
shown  in  Figure  8. 

STDIS 

SHOW  *  ’GAGE(ir  REM1  =  ’LOC:-.U6  m  (PS!)'  DIV  =  6895 
Figure  8.  STDIS  deck  example. 

Note  that  POPT(2)  rr"«t  he  set  equal  to  2  in  the  SINFO  deck  to  activate  the  STDIS  option. 

3.14  SPD1S  Deck.  This  deck  is  similar  to  the  STDIS  deck  except  the  SPDIS  deck  is  for  specifying 
a  desired  output  for  a  parametric  variation  search.  Again,  only  1 1  SPDIS  decks  are  allowed.  In  general, 
this  can  help  identify  which  variables  have  changed  during  a  parametric  variation  and  relevant  related 
output.  The  variables  can  be  selected  from  any  deck  except  SFIND.  SPARA.  SPMAX.  SPDIS  or  STDIS. 
SHOW  identifies  which  variable  is  to  be  printed,  DECK  identifies  the  name  of  the  deck  containing  the 
variable,  and  NTH  corresponds  to  the  number  of  the  deck  (if  there  is  more  than  one).  A  MULT 
or  DIV  of  the  selected  parameter  is  available  to  allow  a  change  in  units.  An  REM1  to  identify  the  column 
of  output  is  also  available.  An  example  of  a  SPDIS  deck  is  contained  in  Figure  9. 
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SPDIS 


SHOW  =  "WEB"  DECK  =  '"PROP"  NTH  =  2 

Figure  9.  SPDIS  deck  example  for  parametric  output  of  the  web  of  the  second  orooellant  deck. 
POPT(5)  must  be  set  equal  to  2,  in  the  5  INFO  deck,  for  this  option  to  be  activated. 

3.15  SFIND  Deck.  The  SFTND  deck  instructs  the  code  to  optimize  a  desired  output  variable  (from 
the  SOUT  list)  through  a  search  over  a  selected  input  parameter.  The  desired  output  variable  is  specified 
by  OUTV  and  the  paremeter  to  be  searched  over  is  defined  by  VARY.  The  location  of  the  VARY 
parameter  is  defined  by  DECK  and  NTH  (if  there  is  more  than  one).  The  initial  value  for  the  parameter 
to  be  varied  is  defined  by  FROM.  Bounds  can  be  placed  on  the  VARY  parameter  with  the  MIN  and 
MAX  variables.  EPS  (incorrectly  defined  in  the  "User  s  Guide")  is  the  smallest  allowable  change  in  the 
parameter  to  be  varied. 

CODE  must  be  set  to  0  to  achieve  a  desired  value  (defined  by  VAL)  of  OUTV.  If  CODE  is  set  equal 
to  1.  the  code  will  attempt  to  maximize  the  value  of  OUTV. 

Up  to  six  SFIND  decks  are  permissible.  If  multiple  SFIND  decks  arc  used,  the  parameter  MULT  is 
a  user-defined  weighting  of  each  rind  deck. 

The  SFIND  deck  has  a  problem  in  selecting  the  correct  increment  to  change  the  VARY  parameter. 
This  results  in  the  code  overshooting  the  exact  point  of  VAL  the  user  desires,  altliough  by  adjustment  of 
EPS.  usable  results  can  be  found.  This  problem  may  be  much  more  critical  if  one  is  using  the  SFIND 
deck  to  maximize  an  output  variable  (OUTV)  with  CODE  equal  to  one. 

3- 16  SPARA  Deck.  A  SPAR  A  search  routine  performs  a  systematic  variation  of  any  chosen  input 
variable,  defined  by  VARY.  The  input  deck  containing  the  variable  to  be  changed  must  be  specified  with 
the  DECK  parameter.  If  there  arc  multiple  decks,  NTH  is  set  to  identify  which  deck  contains  the  variable 
to  be  varied.  The  initial  value  of  the  variable  to  be  tried  is  set  by  PROM,  die  increment  of  change  by  BY. 
and  the  final  value  by  TO. 
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Up  to  four  SPARA  decks  arc  simultaneously  permissible,  which  means  a  four-dimensional  variation 
matrix  can  be  filled.  A  multidimensional  parametric  variation  can  result  in  many  hundreds  or  thousands 
of  runs,  therefore  the  user  must  clearly  identify  the  region  of  interest  as  narrowly  as  possible  to  minimize 
excessive  run  lime. 

The  print  options  for  the  parametric  variations  are  selected  by  P0PT(5)  in  the  5 INFO  deck. 

3.17  $PMAX  Deck.  The  SPMAX  option  performs  an  internal  search  over  a  selected  propellant 
characteristic  to  achieve  a  maximum  breech  pressure.  To  use  this  option,  the  user  must  specify  the 
propellant  variable  to  be  varied  with  the  VARY  parameter  and  give  two  initial  guesses,  TRY1  and  TRY2. 
The  following  guesses  arc  based  on  interpolation  from  these  values.  If  multiple  propellant  decks  are 
present  the  NTH  parameter  must  be  set  to  identify  which  propellant  deck  is  to  undergo  the  changes.  The 
user-defined  maximum  breech  pressure  is  PMAX  and  the  allowed  error  is  EPS.  Bounds  can  be  set  on  the 
range  of  variations  with  the  MIN  and  MAX  variables.  The  LOOP  parameter  sets  a  limit  on  the  maximum 
number  of  iterations  allowed  in  the  search  (default  is  20). 

The  user  can  minimize  convergence  time  of  the  search  by  adjusting  LOOP  and  EPS.  If  there  are 
multiple  local  pressure  maximums.  the  specific  peak  can  be  specified  with  the  NPMX  parameter.  If 
intermediate  information  during  the  search  is  desired,  the  parameter  PRNT  can  be  set  to  1. 

4.  TEST  CASES 

Following  the  conclusions,  four  sample  calculations  are  performed  using  metric  units.  The  first  is  a 
simplification  of  the  standard  test  case  number  6  from  the  IBHVG2  "User’s  Guide."  It  is  simplified  by 
merging  the  two  propellant  decks  into  one,  and  removing  the  multiple  SPMAX  and  SPARA  decks,  which 
iterate  for  a  specific  pressure  and  parametrically  vary  the  propellant  geometry. 

The  resulting  test  deck  is  a  representative  example  of  an  optimized  120-mm  gun  loaded  with  JA2. 
The  first  three  sample  calculations  show: 

1.  performance  of  the  unmodified  gun; 

2.  adding  1  kilojoule  electric  per  gram  of  JA2;  and 

3.  additionally  designing  the  propellant  for  a  flat  pressure  pulse. 


19 


These  three  cases  are  illustrated  in  Figure  10.  A  fourth  case  shows  how  the  trajectory  output  can  be 
printed  in  English  units. 


Figure  10.  Breech  pressure  (MPa)  vs.  time  (ms). 

Case  1  shows  a  typical  propellant  bum.  with  the  rising  pressure  curve,  peaking  at  about  3.5  ms.  Due 
to  the  relatively  long  time  involved  in  this  pressure  rise,  the  projectile  has  time  to  accelerate  down  the 
barrel  before  the  propellant  has  an  opportunity  to  do  its  most  efficient  work.  This  reduces  the  gun 
efficiency. 

In  Case  2.  the  electrical  energy,  which  is  only  16.5%  of  the  total,  is  delivered  in  a  trapezoidal  pulse 
during  the  first  millisecond  of  the  shot,  to  speed  up  the  propellant  combustion.  The  result  is  an  early, 
steep  pressure  rise,  which  increases  the  muzzle  kinetic  energy  by  46%,  but  also  overpressures  the  gun. 

In  Case  3.  the  propellant  gas  generation  rate  is  adjusted  to  restrain  the  maximum  pressure  so  the  gun 
is  not  overpressured.  The  muzzle  kinetic  energy  improvement  of  27%  is  less  than  for  Case  2,  but  still 
well  above  Case  1.  In  this  age  of  "designer''  propellants,  the  gas  generation  rate  versus  time  plot, 
combined  with  the  pressure  history  from  the  simulation,  can  provide  the  guidance  to  the  propellant 
designer  to  construct  a  layered  propellant  that  meets  the  required  performance  specifications. 
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5.  CONCLUSIONS 


IBHVG2  is  capable  of  modeling  a  wide  variety  of  solid  propellants  of  interest  to  the  propulsion 
community  with  relative  ease  of  use.  The  input  file  is  capable  of  specifying  and  activating  all  available 
options  for  a  simulation.  Once  the  propellant  deck  is  satisfactorily  complete,  variations  in  loading  density, 
gun  dimensions,  and  other  optimizations  are  automated  and  relatively  simple  to  make.  With  the 
modifications  to  allow  the  addition  of  parameters  important  in  ET  processes,  modeling  with  1BHVG2  has 
been  extended  to  include  ETC  guns.  It  is  hoped  that  experimental  thrusts  can  be  accentuated  and  directed 
from  results  obtained  using  1BHVG2.  In  addition,  safety  in  the  laboratory  may  be  increased  by  finding 
overpressure  situations  for  a  specific  gun  fixture. 

The  authors  hope  that  this  tutorial  will  be  an  aid  to  the  new  user  of  IBHVG2.  Comments  or 
suggestions  for  a  future  version  are  welcome. 
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APPENDIX  A: 

SUMMARY  OF  INPUT  DECKS 
(Updated  Version  of  Appendix  in  BRL-TR-2829) 
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A.  SYNTAX  RULES 


Input  data  stream  for  IBHVG2  consists  of  a  scries  of  ’’decks."  The  start  of  each  deck  is  identified  by 
a  line  starting  with  a  in  the  first  column  and  followed  by  a  four  letter  mnemonic.  All  lines  of  input 
data  between  decks  should  in  principle  be  closely  related  ballistic  parameters. 

Inputs  are  free  format  in  that  more  than  one  can  be  included  per  line.  Comments  may  be  incorporated 
by  preceding  them  with  a  in  any  column  other  than  the  first  The  formats  are 

variable  rvalue 


for  ordinary  unsubscripted  variables  and 

variable  (sub)  =value,.  -  ,  valuen 

for  subscripted  variables,  with  commas  and/or  blanks  between  consecutive  variable/value  constructs. 
Integer  and  real  formats  are  legal  for  all  numerical  values.  Character  strings  must  be  delimited  with  either 
apostrophes  or  quotes.  All  keywords,  to  include  deck  cards  and  variable  names,  must  be  in  uppercase. 

The  first  blank  between  consecutive  variable/value  fields  is  the  legal  delimiter,  any  additional  ones 
are  ignored.  Blanks  embedded  in  variable  names,  subscripts,  or  numeric  values  are  illegal.  Leading 
blanks  in  character  strings  are  squeezed  out. 

An  omitted  "(sub)"  in  a  subscripted  variable  is  taken  to  be  1.  Consecutive  values  fill  locations  sub, 
sutHl.  sub+2.  etc.,  while  pairs  of  commas  with  nothing  or  blanks  between  them  advance  to  the  next 
subscript  (i.e.,  a  subscript  is  skipped  over).  All  character  variables  default  to  blank  strings  while  numerical 
variables  default  to  either  zero  or  a  convenient  value  listed  below.  If  a  variable  is  multiply-defined,  the 
last  input  overrides  all  previous  ones.  All  syntax  errors  are  flagged.  A  list  of  all  valid  deck  names 
follows  as  well  as  descriptions  of  all  variables  within  each  deck. 

IBHVG2  Version  500  scries  is  a  conversion  to  metric  units  from  the  previous  level  of  IBHVG2 
(Version  408)  which  uses  Imperial  units  of  pounds,  inches,  feet,  and  psi.  In  Version  500  series,  all 
measurements  of  length  are  based  on  the  meter  (m),  all  mass  measurements  on  the  kilogram  (kg),  and 
pressure  input  units  are  the  mcga-pascal  (MPa).  A  list  of  the  constants  used  is  in  subroutine  DATTNP. 
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There  are  a  tew  discrepancies.  Units  of  pressure  in  the  FORTRAN-codcd  equations  use  Pascals  (Pa), 
while  the  unit  of  choice  for  input  is  the  Mega-Pasc.M.  Therefore,  most  instances  of  pressure  during  the 
input  phase  arc  convened  internally  from  MPa  to  Pa  without  the  user’s  knowledge.  Three  instances 
remain  distinct  and  must  be  recognized  by  any  program  user  attempting  a  FIND  process  or  outputs  via 
PDIS  and  TDIS.  In  these  three  cases,  any  reference  to  a  pressure  value  from  the  STRAJ  list  assumes  that 
pressures  will  be  used  in  units  of  Pascals.  PDIS  and  TDIS,  as  special  cases  for  dealing  with  output 
displays,  have  provision  for  multiplication  (MULT)  or  division  (DIV)  of  the  output  values.  The  FIND 
deck  as  an  input  for  an  internal  ceairh  process,  has  no  such  factor  any  reference  to  a  pressure  value  from 
the  STRAJ  deck  must  be  in  Pascals. 

B.  Decks  in  Alphabetical  Order 

SCOMM 

All  cards  between  a  SCOMM  and  the  next  deck  card  are  ignored  to  allow  the  tagging  of  a  data  input 
file  with  user  comments. 

SEND 

This  control  card  signifies  end-of-input  for  the  current  case.  IBHVG2  will  then  execute.  Afterwards, 
it  will  start  reading  the  next  card,  if  any.  This  is  to  allow  the  processing  of  other  runs  or  for  the  code 
itself  to  generate  runs  internally.  An  end-of-input  condition  for  a  user’s  input  deck  is  processed  as  if  the 
SEND  card  was  read  in.  A  frequent  mistake  by  novice  users  is  to  include  the  SEND  and  follow  it  by  a 
blank  line.  The  code  assumes  this  is  another  run  and  proceeds  to  produce  errors  due  to  insufficient 
information. 

$ETC 

Allows  introduction  of  electrical  energy  into  the  combustion  chamber  for  electrothermal-chemical 
(ET( ')  gun  simulation.  This  deck  also  allows  the  gas  generation  rate  to  be  specified  as  a  piecewise  linear 
function  of  time.  This  option  overrides  the  calculation  of  the  surface  areas  and  other  features  of  the 
chemistry. 
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NPWR 

number  of  power-dme  f»'.i  >  (min  =  2.  max  =  20) 

TPWR 

time  of  power  input  (s) 

PWR 

power  at  time  TPWR  (watts) 

NBRN 

numuer  of  gas  generation  rate-time  pairs  (min  =  2,  max  =  20) 

TBRN 

time  of  gas  generation  (s) 

BRN 

gas  generation  rate  at  time  TBRN  (kg/s) 

SFIND 

For  inputs  to  a  general  variation-and-search  al^^dim.  utilizing  function  minimization  techniques. 
Up  to  six  SFIND  decks  may  be  submitted  in  one  run  avowing  a  variation  in  six  dimensions.  A  list  of 
additional  option  variables  for  outputting  is  found  at  the  end  of  this  section. 


VARY 

namr  of  parameter  including  any  subscript 

DECK 

name  of  dtex  containing  parameter,  four  characters  max;  may 
not  be  PARA.  PDIS.  FIND,  PMAX.  or  TDIS 

NTH 

number  of  deck  if  there  are  several  with  same  name 
(default  =  1  ] 

FROM 

initial  value  of  parameter  in  proper  units 

EPS 

convergence  criteria  for  the  varying  parameter 

OUTV 

name  of  desired  output  variable  from  run  completion  variables 
(see  end  of  this  section) 

CODE 

0  to  achieve  desired  value  of  OUTV  variable 

1  to  maximize  OUTV  variable 

VAL 

desired  value  if  CODE  is  0,  else  ignored 

MULT 

multiplier  for  function-minimizer  residual  (default  =  1.0] 

MIN 

minimum  allowable  value  of  VARY  parameter  [default  =  0.0] 

MAX 

maximum  allowable  value  of  VARY  parameter  (default  = 
1.0E+10] 
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SGUN 


For  variables  related  to  gun-tube  geometry,  namely: 


NAME  name  of  gun;  28  characters  max 
TYPE 

CHAM  chamber  volume  [m3] 

CV 

CVOL 

GRVE  groove  diameter  (m) 

LAND  land  diameter  (m] 

TRAV  travel  to  shot-exit  (m) 

LENG 

G/L  ratio  of  groove  to  land  surface  area;  smooth-bore  if  G/L=0 
[default] 

TWST  rilling  twist  [calibers/tum]  [default  =  999] 

LOPT  0  [default]  to  ignore  LDEN 

1  to  calculate  CHAM  from  total  charge  weight  and  LDEN 

2  to  calculate  total  charge  weight  from  CHAM  and  LDEN; 
primer  and  charge  weights  must  be  specified;  they  are  scaled 
proportionally  to  sum  to  required  total  charge  weight 

LDEN  ratio  total  charge  weight/chamber  volume  kgrin3;  used  when 
LOPT  *  1  or  2  [default  *  0.2] 

CAUTION:  LDEN  calculations  are  done  after  C,  ch 
may  have  altered  charge  weights)  is  used. 

CLEN  effective  chamber  length  [in]  for  scaling  in-chamber  p.~ssure 
gauge  locations  [default  «  OTbore  area] 

NGAG  number  of  gauge  locations  (default  =  0,  max  =>  30] 

GLOC  gauge  location  array  of  size  NGAG  [m]  [defaults  =  0.0]. 

Distance  is  measured  (+)  downtube  from  the  breech  or  (-)  into 
chambe  from  the  initial  position  of  the  projectile  base. 
CAUTION:  IBHVC2  will  discard  duplicate  or  out-of-range 
values  and  will  rearrange  distances  in  ascending  order,  if 
necessary.  This  should  be  kept  in  mind  when  using,  for 
example,  a  SSAVE  deck — the  position  you  reference  may  not 
be  the  one  the  computer  code  decided  upon. 
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SHEAT 

For  heat-loss-related  variables,  namely: 


TSHL 

tube  shell  thickness  for  heat  sink  [m]  [default  =  0.0001016] 

CSHL 

shell  specific  heat  (J/kg-K]  covering  a  broad  range  of  steels 
[default  s  460.316318] 

RSHL 

shell  material  density  (kg/m3)  (default  =  7861.0916) 

TWAL 

initial  wall  temperature  [K]  (default  =  293.0] 

HO 

free  convective  heat  transfer  coefficient  for  air  in  the  tube 
[W/m2-K]  (default  =  11.34821852] 

HL 

0  to  ignore  heat  losses  in  energy  balance 

1  [default]  to  include  heat  losses 

SINFO 


For  njn-telated  inputs,  the  variables  are: 


a 

max  integration  time  step  (s]  [default  =  0.00C1] 

n 

1  (default]  for  DELP  in  units  of  time  (s) 

2  for  DELP  in  units  of  projectile  travel  [m] 

DELF 

integrator  logout  and  print  step  (sj  or  |m);  reset  to  DELT  if  DLPU  =  1  and  DELP 
<  DELT 

SOPT 

0  to  suppress  file  storage  of  run  output  [default] 

1  to  write  trajectory  data  for  each  run  into  output  file  STORE  for  post- processing 

2  to  write  only  single-line  summaries  for  each  parametric  run  into  file  STORE; 

ignored  if  nonparametric  runs  are  being  conducted 

GRAD 

1  for  Lagrange  gradient  (default] 

2  for  Pidduck-Kcnt  gradient 

3  Chambrage  gradient.  Requires  description  of  chamber  geometry  in  $GUN  in 
form  of  pairs  of  chamber  diameter  DIAM  versus  distance  from  breech  face 

DIST  with  CPTS  points  in  all 

UNIT 

unit  system  for  output;  not  implemented  yet,  but  will  be  a  choice  between  0  for 

SI  and  1  for  English 
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POPT 

print  option  array  of  size  6  (all  defaults  =  1 ).  Detailed  descriptions  given  below 
between  the  double  bars: 

POPT(l) 

0  to  suppress 

1  to  prim 

POPT(2) 

0  to  suppress  trajectory  print 

1  to  print  default  trajectory  variables 

2  to  print  user-specified  STD  IS  variables 

POPT(3) 

0  to  suppress 

1  to  print  IB  summary 

POPT(4) 

0  to  suppress  blowdown  calculation 

1  tube  blowdown  (to  include  tube  recoil  when  the  recoil  options  is  in  effect,  see 

SRECO)  until  first  rarefaction  wave  reaches  breach  face 

2  tube  blowdown  until  specified  breech  pressure  BLRP  in  SINFO  has  been  reached 
(reduced  printing) 

3  tube  blowdown  as  in  (2)  but  with  more  detailed  intermediate  printing 

NOTE:  Analysis  for  POPT(4)  a  2  or  3  is  an  extension  of  Con.-.r’s  equations  in 

POPT(4)  =  l,  although  this  analysis  is  strictly  valid  only  until  rarefaction  wave  reaches 
breech  face. 

POPT(5) 

0  to  honor  above  print  options  for  every  run  of  parametric  variation 

1  to  honor  above  print  options  for  first  run  of  parametric  variation  and  print  a  single- 
line  summary  thereafter 

2  like  1,  but  summary  print  variables  supplied  via  $PDIS  specifications 

POPT(6) 

is  currently  unused 

RUN 

TTTL 

run  title  on  output  pages:  48  characters  max 

EPS 

maximum  error  for  integrator  time-step  adjustment  and  transition  tolerances  (default  = 
0.0002) 

CONP 

0  for  usual  non-constant  pressure  run  (default] 

1  for  run  with  constant  pressure  maintained  by  varying  the  burning  rate  of  the  single 
charge 

2  for  run  with  constant  pressure  maintained  by  varying  the  surface  area  of  the  single 
charge 

PRES 

desired  constant  brccch  pressure  (psi)  when  CONP  option  is  1  or  2 

BLPR 

desired  brccch  pressure  for  blowdown  calculation  with  POPT(4)  =  2  or  3  in  SINFO 

TOL 

error  tolerance  (MPa)  for  PRES  (default  =  1.0) 
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SPARA 


For  parametric  variations;  up  to  four  SPARA  decks  per  run  permit  a  four-dimensional  matrix  to  be 
systematically  tried.  The  DECK  must  contain  a  nominal  value  of  the  variable  to  be  VARY’d. 


VARY 

name  of  parametric  variable  including  any  subscript 

DECK 

name  of  deck  containing  variable;  four  chan  max;  may  not  be  PARA,  PDIS,  FIND, 
PMAX,  or  TDIS 

NTH 

deck  number  if  there  is  more  than  one  [default *  1] 

FROM 

initial  value  of  variable 

final  value  of  variable 

!■ 

increment/decrement  value 

SPD1S 

Each  deck  defines  one  variable  to  be  printed  in  lieu  of  the  default  set  for  each  line  of  the  parametric 
summary  prim.  There  can  be  up  to  1 1  SPD1S  decks  in  effect:  1  for  each  variable  printed  for  the  interior 
ballistic  cycle.  Be  sure  to  include  the  line  POPT(5)»2  (print  option)  in  the  S1NFO  deck.  To  write  the 
summary  information  to  a  file  attached  to  UNIT=7,  include  the  line  SOPT=2  in  the  SINFO  deck.  (SPDIS 
is  similar  to  the  STDIS  deck.)  A  list  of  additional  option  variables  for  outputting  is  found  at  the  end  of 
this  section. 


SHOW 

name  of  variable  to  print;  four  chars  max 

DECK 

name  of  deck  containing  the  desired  variable;  any  deck  may  be  named  except  PARA. 
PDIS.  FIND,  and  PMAX 

NTH 

number  of  deck  if  more  than  one  (default  =  1] 

MILT 

number  to  multiply  data  value  by  (default  =  1.0] 

DIV 

number  to  divide  data  value  by  (default  =  1.0] 

REMI 

REMK 

20  character  remark  string  to  identify  the  variable  being  displayed 
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SPMAX 


For  variation  of  charge  weight  or  web  to  achieve  a  desired  maximum  breech  pressure.  If  web  is 
varied,  grain  ratios  rather  than  grain  dimensions  may  be  the  better  choices  in  the  SPROP  deck  concerned. 
The  SPROP  deck  must  contain  a  nominal  value  of  charge  weight  and  web.  even  though  one  or  the  other 
will  be  varied  in  the  SPMAX  operation.  Note  that  a  DECK=’PROP*  card  is  not  needed  in  a  SPMAX;  the 
program  knows  that  it  must  be  varying  some  propellant  characteristic.  If  you  include  DECK*' PROP’ . 
the  program  will  complain  but  will  still  do  the  proper  things. 


VARY 

variable  name  in  a  SPROP  deck  including  any  subscript 

NTH 

number  of  SPROP  deck  (default  =  1 1 

TRY1 

first  value  of  VARY  to  try 

TRY2 

second  value  of  VARY  to  try;  third  and  subsequent  guesses  are  based  on 
interpolation.  NOTE:  the  last  two  guesses  from  the  previous  run  are 
employed  as  the  first  two  guesses  in  ixcond  and  subsequent  SPARA  runs 

PM  AX 

maximum  breech  pressure  (MPa]  sought 

EPS 

error  tolerance  (MPa]  for  PM  AX;  (default  *  1.0) 

LOOP 

max  number  of  tries  before  quitting  (default  =  20] 

MIN 

minimum  allowable  value  of  VARY  parameter  (default  *  0.0] 

MAX 

maximum  allowable  value  of  VARY  parameter  (default  =  1.0E+10] 

NPMX 

0  (default|  if  PMAX  refers  to  the  max  breech  pressure  ever  achieved  during 
the  run 

n  where  n  e  ( 1.2.....5)  if  PMAX  refers  to  the  NPMX-th  local  breech  pressure 
maximum  achieved 

SPRIM 


For  primer  data.  The  "primer"  is  considered  completely  burned  at  the  start  of  integration.  For  this 
reason,  it  is  considered  wise  to  include  only  enough  primer  to  reach  sufficient  pressurization  to  ignite  the 
propelling  charge.  Real  primers  arc  usually  modelled  with  1BHVG2  by  simulating  the  function  with  both 
a  SPRIM  and  a  SPROP  deck.  A  typical  correspondence  is  10%  by  weight  for  the  SPRIM  and  90%  as  a 
SPROP.  Constant- pressure  runs  may  include  or  exclude  a  primer.  A  primer  is  mandatory  for  conventional 
simulations. 


34 


NAME 

TYPE 

name  of  primer.  28  characters  max 

GAMA 

specific  heat  ratio 

FORC 

force  (JAgl 

COV 

covolume  (m3/kg] 

TEMP 

flame  temperature  [K] 

CHWT 

WT 

CHGW 

C 

CW 

weight  (kg] 

SPROJ 


For  projectile-related  variables. 


NAME 

TYPE 

projectile  designation;  28  characters  max 

PRWT 

WT 

projectile  weight  [kg] 

COPT 

0  (default]  to  ignore  C/M 

1  to  calculate  PRWT  from  total  charge  weight  and  C/M 

2  to  calculate  total  charge  weight  from  PRWT  and  C/M;  primer  and  all  charge  weights 
must  be  specified,  but  they  will  be  scaled  proportionally  to  sum  to  the  required  total 
charge  weight 

C/M 

used  when  COPT=l  or  2;  ratio  of  total  charge  weight  and  projectile  weight  (default  * 

1.0) 

SOPT 

0  (default]  bypass  PRWT  calculation  based  on  subprojectile  parameters  and  sabot 
formula 

1  find  PRWT  based  on  projectile  weight  estimation  formula  of  Bums 

2  find  PRWT  based  on  same  weight  estimation  formula  using  coefficients  yielding  15% 
lighter  sabot 

WTSP 

subprojectile  weight  (kg] 

USP 

subprojectile  length  (m) 

DSP 

subprojectile  diameter  (m] 

PD&S 

max  design  pressure  (MPa)  for  saboted  projectile 

SABO 

sabot  weight  (kg];  not  input— set  by  IBHVG2 
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SPROP 


Defines  a  main  propelling  charge  clement  16HVG2  will  recognize  up  to  five  such  decks  and  considers 
them  independently  (i.e.,  order  is  unimportant).  The  following  are  basic  input  variables  most  useful  for 
describing  homogeneous,  undeterred  grains: 


NAME 

TYPE 

name  of  propellant:  28  chan  max 

RHO 

DENS 

density  [kg/m3] 

GAMA 

specific  heat  ratio 

FORC 

force  [J/kg] 

COV 

covolume  (m3/kg] 

TEMP 

flame  temperature  [K] 

CHWT 

WT 

CHGW 

C 

CW 

weight  (kg] 

ALPHA 

burning  rate  exponent,  a 

BETA 

burning  rate  coefficient,  p  (m/s  -  MPa®]  (where  burning  rate  «pp®  (m/s)  where  P  is 
mean  pressure  in  (MPa] 

EROS 

erosive  burning  coefficient,  empirical  factor  multiplied  by  projectile  velocity  [m/s]  to 
add  to  burning  rate  (default  =  0.0] 

GRAN 

FORM 

granulation  code  chosen  from  7PF  (or  7P),  1PF  (IP),  CORD,  RECT  (SLAB),  SPHR 
(BALL),  SLOT,  37HX  (37H),  19HX  (19H),  19PF  (19P),  GEN,  PIE  (STAR).  GHEX 
(HEX).  MONO 

The  following  form  functions  have  been  added  to  IBHVG2  since  the  User's  Guide  (Anderson  and 
Fickic  1987)  was  published. 

ai  Rattened  ball  powder  grains  produced  from  ball  powder  by  rolling  to  a  pancake  shape  with 
constant  thickness.  The  form  function  name  is  PAN  or  CAKE,  and  the  geometry  is  defined  by 
DIAM  and  THCK.  DIAM  corresponds  to  (he  outer  diameter  or  maximum  width.  THCK 
corresponds  to  the  thickness  of  the  rolled  portion  or  distance  between  parallel  plane  surfaces.  If 
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only  (he  original  ball  diameter  and  target  thickness  arc  known,  the  fotm  function  will  calculate  the 
final  disk  diameter  when  DIAM  is  the  negative  of  the  value  of  the  original  ball  diameter  (DIAM 
<  0).  In  this  ease,  the  computed  disk  diameter  is  printed  during  the  propellant  description  phase 
of  program  output. 

b)  Rosette  form  function  ROSE.  The  geometry  is  similar  to  that  described  in  ARBRL-MR-03380. 
The  required  dimensions  are:  perf  diameter,  web  thickness  (as  WEB  or  as  W|  and  WO  where  WO 
corresponds  only  to  outermost  web),  and  length.  Additionally,  the  number  of  perf  rings  NRNG  as 
described  in  the  general  hexagon  form  function  must  be  specified.  Layering  is  allowed  on  any  or 
all  surfaces. 

c)  General  hexagon  (GHEX)  form  function.  GHEX  may  have  a  different  thickness  of  the  outer  web. 
Previously,  all  webs  were  equal. 


FORM 

PAN(CAKE) 

ROSE 
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wi 

WIN 

WO 

WOUT 

WM 

WMID 

WEB 

WB 

D 

DIAM 

GDIA 

PD 

DP 

PDIA 

SLOT 

NSLT 

NRNG 


L 

GL 

LEN 

GLEN 

WDTH 

THCK 

NSUR 

DEPB 

SURF 

IGNC 


inner  web  for  cylindrical/hcxagonal  grains  [ml 

outer  web  for  cylindrical/hexagonal  grains  [m] 

middle  web  for  cylindrical/hexagonal  grains  (mj 

common  value  of  all  inner,  middle,  and  outer  webs  (m  j;  also  resets  WI/O 

grain  diameter,  if  applicable  (mj 


perf  diameter,  if  applicable  [m] 
(default  =  0] 


slot  width  in  SLOT  and  PIE  grains  [m] 
number  of  diametral  slots  in  PIE  grain  2  2 

number  of  concentric  rings  of  perfs  around  central  perf  in  HEX  and  rosette  grains;  2 0 
[default  =  2] 

grain  length,  if  applicable  [m] 


grain  width  in  RECT  grains  (mj 
grain  thickness  in  RECT  grains  (m) 

number  of  depth/surface  pains  for  GEN  grain  [default  =  1.  max  =  20] 

depths-bumed  array  (m)  for  GEN  grains;  ignored  if  NSUR  *  1 

surface  area  array  [m2]  for  GEN  grains 

code  specifying  charge  ignition: 

0  to  ignite  at  start  of  integration  (default) 

1  to  ignite  at  some  time  (sj 

2  to  ignite  at  some  projectile  travel  (m] 

3  to  ignite  at  some  mean  pressure  (MPa) 

4  to  ignite  at  some  Z  (mass  fraction  burned)  of  the  charge  described  by  the  previous 
SPROP  deck 
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THRC 

threshold  value  for  ignition  if  IGNC  >  0 

IGNS 

array  of  size  3  specifying  ignition  codes  for  the  perf,  end,  and  lateral  grain  surfaces, 
respectively;  same  codes  as  IGNC;  ignored  if  IGNC  >  0.0 

THRS 

array  of  threshold  values  for  surface  ignition 

PA-B 

0  (default]  for  standard  treatment  of  charge;  set  automatically  if  grain  has  no  perfs  or  if 
runs  has  either  constant-pressure  option  set 

1  for  perf-augmented  burning  of  charge  until  grain  fracture 

DSCF 

Robbins-Horst  discharge  coefficient  in  weight  flux  computation  for  perf-augmented 
burning  model  [default  =  1.0] 

FRAC 

0  [default]  for  no  user-defined  grain  fracture  criterion  of  charge  with  PA-B  =  1.  NOTE: 

web  bum-through  turns  off  perf-augmented  burning. 

1  to  specify  pressure-difference  grain  fracture  threshold  for  charge  with  PA-B  =  1. 

NOTE:  web  bum-through  triggers  fracture. 

NFRG 

number  of  fragments  when  1PF  grain  fracture  2  1 

THRF 

absolute  value  of  difference  between  chamber  and  perf  pressures  for  grain  fracture 
[MPa];  ignored  if  FRAC  =  0  (default  =  0.0] 

IBHVG2  will  calculate  grain  dimensions  if  given  either  actual  measurements,  or  alternatively,  a  single 
measurement  combined  with  ratio  specifications.  This  latter  technique  is  quite  useful  when  the  code  is 
used  to  optimize  a  charge  design.  For  example,  a  CORD  grain  can  be  defined  by  D  and  L,  or  by  L  and 
the  ratio  L/D.  For  a  7-perforated  grain.  WI  and  DP  toother  with  WI/O,  D/DP,  and  L/D  "/ill  completely 
specify  the  grain  geometry.  The  slot  width  of  a  SLOT  or  PIE  grain  will  be  calculated  given  SW/D  and 
D.  Care  must  be  taken  to  avoid  over-specified  or  inconsistent  information  since  the  computer  program 
may  make  arbitrary  or  poor  decisions.  1BHVG2  prints  the  grain  dimensions  prior  to  ballistic  calculations; 
experienced  users  of  the  code  always  scrutinize  this  portion  of  the  output  for  any  surprises. 


n 

grain  length/grain  diameter 

grain  length/perf  diameter 

m 

grain  diamctcr/pcrf  diameter 
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L/WD 

grain  length/grain  width 

WD/T 

grain  width/grain  thickness 

SW/D 

slot  width/grain  diameter 

WI/O 

inner  web/outer  web  for  multiperforated  grains;  web  eccentricity 
for  single  perf  grains  (default  a  1.0] 

For  deterred  grains,  further  (or  replacement)  inputs  to  handle  variable  thermochemistry  and  burning 
rates  are  necessary.  Visualize  each  grain  having,  in  general,  three  disjoint  surfaces: 


Surface 

Description 

P 

E 

L 

perf  (includes  the  slot,  if  any) 

end 

lateral 

Extending  inward  from  each  surface  are  layers  1.  2,  3,  and  4  (the  innermost  or  "core"),  and  any  or 
all  of  (he  first  three  may  be  of  zero  thickness  (nonexistent)  on  one  or  more  of  the  P.  E,  or  L  surfaces. 
The  core  layer  is  always  present  and  its  properties,  unlike  those  of  the  outer  layets,  are  independent  of 
depth.  As  there  is  only  one  common  4th  layer,  L-specifications  for  the  core  override  all  others  to  prevent 
conflicts.  All  depth-varying  properties  are  consecutively  defined  at  the  outer  surface  of  each  layer  from 
the  outside  in  so  each  of  the  following  is  an  array  of  size  4.  Linear  interpolation  is  performed  by 
IBHVG2  when  intermediate  values  are  required. 


RHOP 

RHOE 

RHOL 

density  (kg/m3] 

on  P  surfaces 
on  E  surfaces 
on  L  surfaces 

GAMP 

GAME 

GAML 

specific  heat  ratio 

on  P  surfaces 
on  E  surfaces 
on  L  surfaces 

FRCP 

FRCE 

FRCL 

force  ( J/kg] 

on  P  surfaces 
on  E  surfaces 
on  L  surfaces 

COVP 

COVE 

COVL 

covolumc  (m3/kg] 

on  P  surfaces 
on  E  surfaces 
on  L  surfaces 

TMPP 

TMPE 

TMPL 

flame  temperature  [K] 

on  P  surfaces 
on  E  surfaces 
on  L  surfaces 
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Array  of  size  3  specify  the  transition  depths  (m)  between  the  first  and  second,  the  second  and  third, 
and  the  third  and  fourth  layers,  respectively,  for  each  surface  (measuring  from  the  surface). 


Transition  depths  will  be  calculated  by  IBHVG2  if  positive  ratios  specifying  transition  depths  as 
fractions  of  surface-to-surface  distance  are  input.  Arrays  of  size  3  store  the  ratios  [the  default  for  all 
values  is  0.0]. 


DP/S 

ratios  of  transition  depths  on  P  surfaces  to  Wl;  used  for  perforated  grains  only 

DE/S 

ratios  of  transition  depths  on  £  surfaces  to: 

•  min(GL.WDTH)  for  RECT 

•  GL  for  perforated  and  CORD  grains 

DUS 

ratios  of  transition  depths  on  L  surfaces  to: 

•  D  for  CORD  and  SPHR  grains 

•  THCK  for  RECT  grains 

•  WI  for  perforated  grains 

FP/L 

ratios  of  DP/S  to  DUS;  DP/S  values  will  be  found  given  FP/L  and  DL/S;  useful 
for  defining  a  relation  between  transition  depths  on  P  and  L  surfaces 

DEPP 

transition  depths 

on  P  surfaces 

DEPE 

on  E  surfaces 

DEPL 

(defaults  *  0.0] 

on  L  surfaces 

Finally,  burning  rate  for  the  outer  surface  of  each  layer  may  be  incorporated  via  several  alternative 
methods.  It  should  be  noted  that  interpolation  is  linear  for  depths  and  betas,  but  logarithmic  for  alphas 
and  tables  of  pressure  versus  burning  rates. 


NTBL 

method  of  specifying  burning  rate  inputs  ...  absolute  value  is  number  of  table  entries  on 
every  surface  layer,  range  of  values:  -10  to  +10 

<  0  to  specify  betas  and  alphas  as  tabular  functions  of  mean  pressure  (MPa),  NTBL 
triples  in  all 

=  0  to  define  one  bcta/alpha  pair  [default] 

>  0  to  specify  burning  rates  (m/s]  as  tabular  functions  of  mean  pressure  [MPa],  NTBL 
pairs  in  ail.  NOTE:  if  NTBL  =  1,  burning  rate  is  constant  so  corresponding 
pressure  value,  if  any.  is  ignored. 

PRIL 

pressures  [MPa]  on  outside  of  layer  1,  L  surface. 

PR2L,  PR3L,  PR4L,  PR  IE.  PR2E,  PR3E,  PR4E,  PR  IP.  PR2P,  PR3P,  PR4P  are  defined 
similarly. 

BRIL 

burning  rates  (m/s]  on  outside  of  layer  1,  L  surface.  BR2L,  BR3L.  BR4L,  BRIE, 

BR2E,  BR3E,  BR4E.  BR1P,  BR2P,  BR3P,  BR4P  are  defined  similarly. 

CF1L 

burning  rate  coefficients  (betas)  on  outside  of  layer  1,  L  surface.  CF2L,  CF3L,  CF4L, 

CF1P,  CF2P.  CF3P.  CF4P  are  defined  similarly 

EX1L 

burning  rate  exponents  (alphas)  on  outside  of  layer  1 ,  L  surface.  EX2L,  EX3L,  EX4L, 

EX  IE,  EX2E,  EX3E,  EX4E,  EX  IP.  EX2P,  EX3P,  EX4P  are  defined  similarly 

SRcCO 


For  recoil  data.  In  the  model  currently  available,  the  gun  tube  freely  recoils  under  the  influence  of 
breech  pressure  less  resistance  pressure.  This  option  has  rot  been  fully  tested  and  should  not  be 
considered  reliable. 


RECO 

0  for  fixed  tube,  no  recoil  [default] 

1  to  employ  recoil  option 

RCWT 

WT 

weight  (kg]  of  tube  and  recoiling  parts 

NAME 

TYPE 

name  of  recoil  system;  28  diameters  max 

SRESI 


This  precedes  resistance  pressure  inputs.  The  variables  are; 


NPTS 

number  of  travel/prcssurc  pairs  (min  =  0,  max  =  20] 

TRAV 

projectile  travel  army  of  size  20  (m) 

PRES 

resistance  pressure  array  of  size  20  (MPa) 

AIR 

0  to  suppress  adding  in  air  resistance 

1  to  include  air  resistance  (default) 

HTFR 

fraction  of  work  done  to  overcome  barrel  rcsis’ance  which  is  used 
to  prc-hcat  t  be  wall;  0.0  2  HTFR  £  1.0  (default  =  0.0] 
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A  factor  to  allow  th.  be:;  r  ;istance  profile  to  be  raised  or  lowered  during  an  iterative  searching 
process  can  be  specified.  Tnc  factor  works  in  conjunction  with  SFIND  or  SPARA. 


FACT  factor  by  which  the  resistances  (in  MPa)  will  be  multiplied  [default  =  1.0] 

RFPT  point  in  the  PRES  ana,  vwiere  factoring  will  begin:  factor  will  be  applied 

to  all  subsequent  points  [default  =  3] 


JSAVE 

Actually,  not  a  deck  card,  but  a  control  card.  If  it  is  the  first  card  in  a  run  after  an  SEND  card,  the 
values  of  all  input  variables  are  retained,  so  that  succeeding  decks  need  only  update  selected  variables. 
If  an  SEND  card  is  not  followed  by  a  SSAVE  card,  all  input  variables  must  be  reinitialized. 

STDiS 

Each  deck  defines  one  variable  to  be  printed  in  lieu  of  the  default  set  for  each  line  of  the  trajectory 
print.  There  can  be  up  to  11  STDIS  decks  in  effect— one  for  each  variable  printed  for  the  interior  ballistic 
cycle.  Be  sure  to  include  the  line  POPT(2>=2  (print  option)  in  the  SINFO  deck.  To  write  the  trajectory 
information  to  a  file  attached  to  UNIT=7,  include  the  line  SOPT=2  in  the  SINFO  deck.  (STDIS  is  similar 
to  the  SPDIS  deck).  The  list  of  option  variables  available  for  use  is  found  at  the  end  of  this  section. 


SHOW 

name  of  variable  to  print  from  deck  TRAJ;  four  characters  max 

MULT 

number  to  multiply  data  value  by  [default  =  i.0] 

DIV 

number  to  divide  data  value  by  [default  *  1.0] 

REMI 

20  character  remark  siring  to  identify  the  variable  being 

REMK 

displayed 

Trajectory  Variables 

Below  is  a  list  of  keyword  variables  which  can  be  used  in  conjunction  with  the  STDIS,  SFIND,  and 
SPDIS  decks  to  reference  quantities  other  than  the  default  scl  To  use  them  with  SPDIS  and  SFIND,  one 
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must  include  the  line  DECK=TRAJ';  a  reference  within  STDIS  assumes  the  keyword  will  come  from  the 
following  list.  Test  cases  3  demonstrates  the  use  of  STDIS  to  change  the  trajectory  display  printed  during 
the  ballistic  cycle. 


MEAN 

mean  gas  pressure  in  chamber  [Pa] 

PRFP 

array  of  size  5  of  mean  gas  pressure  in  perfs  [Pa]  of  each  charge;  equal  to 
mean  if  no  perf-augmented  burning  in  that  charge  at  current  time  step 

GAGE 

array  of  size  30  of  gauge  pressure  [Pa] 

BRCH 

breech  pressure  [Pa] 

BASE 

pressure  at  projectile  base  [Pa] 

PDOT 

d(mean  chamber  pressure )/dt  [Pa/s] 

PDTP 

array  of  size  5  of  d(mean  perf  pressure )/dt  [Pa/s],  one  for  each  charge 

TBAR 

mean  gas  temperature  in  chamber  [K] 

PRFT 

array  of  size  5  of  mean  gas  temperature  in  perfs  [K]  of  each  charge;  equal 
to  TBAR  if  no  perf-augmented  burning  in  that  charge  at  current  time  step 

FRCR 

bore-friction  resistance  pressure  [Pa] 

AIRR 

air  resistance  pressure  [Pa] 

TOTR 

sum  of  FRCR  and  AIRR 

TWAL 

temperature  of  tube  wall  shell  [K] 

WTB 

array  of  size  5  of  weight  burned  of  each  charge  [kg];  reference  by 
subscript 

WTBR 

array  of  size  5  of  weight-burning  rate  of  each  charge  [kg/s];  reference  by 
subscript 

WTBT 

total  weight  of  gas  in  chamber  [kg] 

PREJ 

projectile  translational  kinetic  energy  [J] 

PRJ% 

PRJE  *  100  /  TOTE 

PRPE 

propellant  and  gas  kinetic  energy  [J] 

PRP% 

PRPE  *  100  /  TOTE 

ROTE 

projectile  rotational  kinetic  energy  [J] 

ROT% 

ROTE  *  100 /TOTE 
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FRTE 

barrel-frictionr.l  work  to  tube  (J| 

FRT% 

FRCE  *  100  /  TOTE 

FREE 

barrel-frictional  work  not  absorbed  as  heat  to  the  tube  wall  (J) 

FRE% 

FREE  *  100  /  TOTE 

DRGE 

work  done  against  air  in  barrel  (J] 

DRG% 

DRGE*  100 /TOTE 

RECE 

kinetic  energy  of  recoiling  tube  [J] 

REC% 

RECE  •  100 /TOTE 

HETE 

energy  lost  as  heat  convected  to  tube  wall  [J] 

HET% 

HETE  *  100  /  TOTE 

LOSE 

sum  of  all  energy  losses  (J] 

LOS% 

LOSE  *  100  /  TOTE 

TOTE 

total  chemical  energy  released  by  combustion  (J) 

EDOT 

d(TOTE)/dt  (J/s) 

GASE 

internal  energy  of  gas  (J),  i.e.,  TOTE  -  LOSE 

GAS% 

GASE  *  100  /  TOTE 

SRF 

array  of  size  5  of  burning-surface  area  of  each  charge  (m3);  reference  by 
subscript 

SRFT 

total  surface  area  of  all  ignited  charges  (m3) 

TIME 

time  [s] 

TRAV 

projectile  displacement  from  initial  position  [m] 

VEL 

ground-based  projectile  velocity  (m/s] 

ACCL 

projectile  acceleration  (G's] 

Z 

array  of  size  5  of  charge  weight  fractions  burned;  reference  by  subscript 

DB-P 

array  of  size  5  of  depth  burned  into  P  surface  of  each  charge  (m); 
reference  by  subscript 

DB-E 

same  as  above  for  E  surfaces 

DB-L 

same  as  above  for  L  surfaces 

DBF 

same  as  above  for  F  (fracture)  surfaces 

BR-P 

array  of  size  5  of  burning  rate  on  P  surface  of  each  charge  (m/s);  reference 
by  subscript 
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BR-E 

same  as  above  for  E  surfaces 

BR-L 

same  as  above  for  L  surfaces 

BR-F 

same  as  above  for  F  (fracture)  surfaces 

EPWR 

electrical  power  (waits) 

electrical  energy  integrated  from  t*0  [J] 

XBRN 

gas  generation  rate  [kg/s] 

RUN  COMPLETION  VARIABLES 

These  are  output  variables  defined  after  a  complete  1BHVG2  run.  Typically  they  represent  some 
global  extrema  which  can  only  be  ascertained  at  the  conclusion  of  ballistic  computation.  For  example, 
a  maximum  pressure  is  determined  after  shot  ejection,  while  there  could  be  several  local  maxima  during 
the  pressure  history.  The  variables  may  be  referenced  by  name  using  $PDIS  (with  DECK=’OUT')  or  by 
SFIND  as  the  value  inserted  for  OUTV 


PMAX 

max  breech  pressure  [MPa] 

HUMP 

array  of  size  5  of  local  breech  pressure  maxima  [MPa] 

GMAX 

array  of  size  30  of  gauge  pressure  maxima  [MPa] 

VMUZ 

muzzle  velocity  [m/s] 

AMAX 

max  acceleration  [g's] 

BMAX 

max  base  pressure  [MPa] 

X@BO 

array  of  size  5  of  projectile  position  [m]  at  charge  burnouts:  reference  by 
subscript 

PMUZ 

base  pressure  at  shot  exit  [MPa] 

ZMUZ 

array  of  size  5  of  charge  weight  fractions  burned  at  shot  exit;  reference  by 
subscript 

IMPL 

final  impulse  (momentum  (N-s|) 

LDEN 

calculated  loading  density  (kg/m3) 
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APPENDIX  B: 
SAMPLE  CASES 
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Intentionally  left  blank. 
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B.l.  SAMPLE  CASE  1 


This  is  a  simplified  version  of  IBHVG2  Benchmark  test  case  6.  The  simplifications  are  the  removal 
of  the  SPMAX  and  SPARA  decks,  and  combination  of  the  two  propellant  decks  into  a  single,  8-kg 
load. 

While  the  input  is  basically  the  same,  the  output  is  not.  The  JTDIS  decks  direct  the  display  of  three 
new  variables:  the  electric  power,  the  integrated  electric  energy,  and  the  propellant  gas  generation  rate. 
The  first  two  of  these  are  0  in  this  case.  The  third  shows  the  gas  generation  rate  as  a  function  of  time. 

This  case  is  a  baseline  against  which  the  other  two  cases  should  be  compared.  The  breech  pressure 
traces  for  all  three  cases  are  shown  in  Figure  10  on  page  20. 


The  input  deck: 


SC04S1 

SIMPLIFIED  I8HVC2  SCNCKMRAK  IUI  CMC  6 


siwro 

MUM  -  '120m  T6  WITHOUT  PRESSURE  CONSTRAINT'  OUT  -  5C-6  OCLP  •  St-5 
CHAD  -  2  PORT  -  1.2,  1,0,2  SOPT  •  0 

CPS  •  0.001 

(HEAT 

T3HL  •  0.000110  CJHL  *  460.310106  AJML  •  7(61.0*16 

TWAL  -  293  HO  -  1 1 . 3482 IS  HI  -  l 

I  GUN 

M M4E  -  '120HH  GUN  TCST  CASC'  CHAM  -  0.00**46*48  CRVE  -  0.11»*S*6 
LAND  -  0.11***»6  S/l  -  1.  TRAV  -  4.7525*4 

7NST  -  »* 

SPROJ 

NAME  -  'APrSUS'  PH NT  -  7,0*072 

scorn 

'POIS'  VALUES  USED  WITH  * AAAMETAJC  PRINT  OPTION  P0PTIS1-2 

spots 

SHOW-'PMAX'  DECK-' OUT'  OJV-6S*6.757 

SPOIS 

SHOW-'CHNT'  DECK-' PROP'  NTH-2  OJV-O. 4535*237 

SPOIS 

SHOW-' 01AM'  OECK-' PROP'  NTH-2  OIV-O.02S4 
SPOIS 

SHON-'PO'  DECK-' PROP'  N7H-2  DIV-0.0254 

SPOIS 

SHON-'MEB'  OECK-' PROP'  HTH-2  DIV-0.0254 
SPOIS 

SHOW-' VMUZ'  OECK-' OUT'  O7V-O.3O40 
SPOIS 

SHOW-' ZMUZ ( 2 ) '  DECK-'  OUT' 

SPOIS 

SHOW-' IDEM'  DECK-' OUT' 
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suit 


HPT 3  -  4  AIR  -  1 

TttAV  -  0,  0.02032.  0.0742,  4.7496 

PRES  •  0.4194757,  17.234S925,  3.4694757,  0.4694757 


STD1S 

SHOM-'  TIME' 

trots 

SKON-'TRAV' 

0tV'.0254  REMK-' INCHES' 

4TDI* 

SHOW-'VEL' 

HULT-3. 2606111  RDW-'Pr/S' 

4TDIS 

SHOW*' ACCL' 

REMR''  GRAVITIES' 

trots 

SHOtt*' 6MH' 

0tV'4694.757  RIM-'PSI' 

trots 

SROW-'MEAM' 

DtV'4694 .757  RDfll-'PSI' 

trots 

SHOW-'tASC' 

OIV'4694.757  REMf'PSI' 

trots 

SHOW-’EPWR' 

REME"*ATTS  ELECTRIC' 

trots 

SMOW-'EIME' 

REMR"  JOOUS  ELECTRIC' 

trots 

tHOW-'t(l)' 

trots 

SHOW''  XSRM' 

REMR-'SURM  RATE' 

tRCCO 

NAME  •  'NONE' 

RE  CO  -  0  RCCT 

tPRtM 

NAME  •  'SEMITE 

’  CHCT  -  0.301571944 

OAMA  •  1.25 

TORC  -  415174.7375 

COV  -  0.001061619  TEMP  •  2000 

tPROP 

kmc  -  '  ja2  7p>  cmrr  -  6.0  cram  -  ' 7pp' 

RHO  -  1564.411646  OAMA  -  1.2257  TORC  -  1142277.932 

COV  -  0. 000992776  TO*  -  1400  CMS  -  0.0000000 

CTRL'S 

PR4L-  13.769514,  27.579026.  46.94757,  172.146925 

•R4L-  0.02447,  0.016406.  0.074422.  0.144424 

LCM  -  0.0141122  0IAM  •  0.010446  PO  -  0.000506 

MI  -  0.0019104  HO  •  0.0016794 


SEND 


Produced  the 


following  output: 
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T<PC:  CUN  TEST  CASE  CHAMBER  VOLUME  (Ml):  0.0099S  TRAVEL  (M)  J  9.7S2S9 
GROOVE  OlAMf  lER  (M)  :  0.11999  LAND  DIAMETER  (M)  :  0.11999  GROOVE/LAND  RATIO  (->1  1.000 
TWIST  iCALS/TURN) t  99.0  BORE  AREA  (M2):  0.01131  HEAT-LOSS  OPTION:  1 
SHELL  THICKNESS  (HI:  O.POOIM  SHELL  CP  (J/KG-K) :  460. 3163  SHELL  DENS  ITT  (KG/M3)  :  7861.0920 


1201*1  16  WITHOUT  PRESSURE  CONSTRAINT  IBHVG2.S04  DATE  TIME 
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!^f*1  T t»  WITHOUT  PHL-iSURE  CONSTRAINT  ISi!VC2.S0«  OATE  TINE 
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LOCAL  PRESSURE  MAX  DETECTED 
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TRAJLCTOftY  VARIABLES 
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12GNW  rt>  WITHOUT  PRESSURE  CONSTRAINT  1RHVCJ.504  DATE  TIME 


B.2.  SAMPLE  CASE  2 


This  case  differs  from  sample  case  1  by  the  addition  of  1  kJ  of  electric  energy  per  gram  of  JA2 
propellant.  For  8  kg  of  JA2,  that  means  8  MJ  of  electrical  energy. 

This  energy  is  (arbitrarily)  delivered  by  means  of  a  trapezoidal  pulse  in  approximately  1  ms.  The 
precise  details  are  shown  in  the  SETC  deck  below.  The  trapezoid  rises  linearly  in  the  first  100  pa,  remains 
flat  for  900  ps,  and  then  drops  linearly  in  a  final  100  ps.  The  two  triangular  rising  and  falling  segments 
may  be  combined  geometrically  to  show  that  the  net  result  is  a  1-ms  pulse. 


This  early  delivery  of  energy  causes  the  JA2  to  reach  higher  pressures  than  it  was  originally  designed 
for,  increasing  the  bum  rate  and  therefore  substantially  increasing  the  peak  pressure  and  decreasing  the 
time  required  to  reach  it. 

The  input  deck: 


team 

SIMPLIFIED  IBHVC2  BENCHMARK  TEST  CASE  « 

WITH  A  1  MS  rLAT  ELECTRIC  PULSE 

CONTAINING  A  TOTAL  OF  1KJ  ELECTRIC  PER  GRAM  Of  PROPELLANT 

nwro 

RUN  -  '  120m  T»  WITH  1  KJ/C  ELECTRIC'  CELT  -  5C-5  DELP  -  5E-5 
GRAD  -  2  PORT  -  1.2. 1,0.2  SOFT  •  0 

CPS  -  0.0S 

SnCAT 

T3HL  •  0.000110  CSHL  •  460.3163186  RSHL  •  7861.0916 

TWAL  -  293  MO  -  11.348218  ML  -  1 

6  GUN 

NAME  •  '120rn  CUN  TEST  CASE'  CHAM  •  0.009946948  CRVE  -  3.1199896 
LAND  -  0.1199896  G/L  -  1.  TRAV  -  4.752594 

TWST  -  99 

9PROJ 

NAME  •  'APrsOS'  PANT  •  7.09872 

team 

'P0I8'  VALUES  USED  WITH  PARAMETRIC  PRINT  OPTION  P0PT(5)-2 
SPDIS 

SHOW-'PMAX'  DECK-' OUT'  OIV-6894.757 
SPDIS 

SHOW-'CHWT'  DECK-'  PROP'  NTH -2  OIV-O. 45359237 
4POIS 

SHOW-'  DIAM'  DECK-' PROP'  NTH-2  0IV-0.02S4 
SPDIS 

SHOW-'PO'  OECK-' PROP'  NTH-2  OIV-0.02S4 
SPDIS 

SHOW-' WEB'  DECK-' PROP'  NTH-2  OIV-0.02S4 
SP01S 

SHOW-' VMUZ'  OECK-'OUT'  OIV-0.3048 
SPDIS 

SHOW-' ZMU2 (2) '  OECK-'OUT' 

SPOIS 

SHOW-'LOEN'  OECK-'OUT' 

9RESI 

NPTS  -  4  AIR  •  1 

TRAV  -  0,  0.02032,  0.0762,  4.7499 

PRES  -  0.6094757.  17.2368925.  0.6894757,  0.6094757 


STDIS 

SHOW-' TIME' 

STDIS 

SHOW-' TRAV'  DIV-.0254  REMK-' INCHES' 
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sxors 


STDIS 

STDIS 

STDIS 

STDIS 

mis 

STDIS 

STDIS 

STDIS 


SHOW-' VEL' 
SHOW-' ACCL' 
SHOW-' BRCH' 
SHOW-'  MEAN' 
SHOW-'  BASE' 
SHOW-'  EPWR' 
SHOW-'  tort' 
SHOW-'  1(1)' 
SHOW-'  JCBRN' 


MULT-3. 2808331  RD4K-' fT/S' 
REMK-*  GRAVITIES' 

01V-8894 .757  REMX-'PSI' 
OIV-S8B4 .757  AEMX-'PSI' 
DIV-SBB4 .757  REMX-'PSI' 
REMR-'  WATTS  ELECTRIC' 
ROW!-' JOULES  ELECTRIC' 

REMX-'SURN  RATE' 


SRCCO 

NAME  •  'NONE' 


RECO  -  0  RCWT  -  0 


SPRIM 


STROP 


NAME  •  '  BENITE'  CHWT  -  0.0015739** 
GAMA  •  1.25  rORC  •  S35176.7375 
COV  •  0 .001083819  TEMP  -  2000 


NAME  -  '  JA2  7p'  CHWT  -  8.0  GRAN  -  '  7pr< 

RHO  •  158*. *11868  GAMA  •  1.22S~  TORC  -  1142277.912 

COV  -  0.000992778  TD4P  -  3400  EROS  -  0.0000000 

NT8L-4 

PR4L-  13.789514.  27.579028.  <8.94757.  172.3*8925 

BR4L-  0. 02**7.  0.038*08.  0.074422.  0.1MC24 

LEN  •  0.01*3322  DIAM  •  0.010**8  70  •  0.000508 

WI  -  0.0019304  W0  -  0.001879* 


SC0M4 

OEriNC  A  TRAPEZOIDAL  ELECTRIC  POWER  PULSE 
S  GIGAWATTS  FOR  1  MILLISECOND. 

WITH  100  MICROSECOND  LINEAR  RAMPS  UP  AND  DOWN. 

JETC 

NPMR  •  4 

PWR  -0.0.  SE9,  8E9.  0.0 

TPWR-O. 000. 0.0001. 0.001. 0.0011 

SCND 


Produced  die  following  output: 
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CARD  1  —  >  SCOW 

CARD  2  —  >  SIMPLIFIED  IBHVG2  BENCHMARK  TEST  CASE  6 

CARD  3  —  >  WITH  A  1  MS  FLAT  ELECTRIC  PULSE 

CARO  4  -->  CONTAINING  A  TOTAL  OF  IKJ  ELECTRIC  PER  GRAM  OF  PROPELLANT 


“fc  5! 

R-O  WO 

a.  ox 


u  *  «  o  o  O  ♦  i 

£  ^  Q  O 

<  •  «"> 

B  •  N  M  *  C 

<  <•  I  •  W  o  * 

Ai  I  X  Z  I  I  < 

>  t  c  •  z  > 

x~zzxh~* 

S  °  .  .  5  5  °  B 

S  *  ft  Ot  *  *  O  * 

h  o  o  •  an*  t* 

q  s  s  c  a  o  s  i  b 

woa.a<Q£oxO 

w  »  •  *  £  CL  •  u  % 

3  i  •  •  a.  *  «  u  i 

XX  If  -  I  X  O  X 

*i  u  u  u  •  x  u  u 

I  1  H  «  j  I  8  S  8  * 

»  S  S  *5  7-358  . 

»  n  s  s  n  5  9  ! 

w*«»»**»*  7 

;A.v)|n|M)in;v);M|«i|vi|  *- 

:%  hQmOwAmQhOwOwOhO  W 

f  a  x  a  x  a  *  a  z  o  x  o  ±  ©  z  o  x  h 

»  a*«iOiVia.wa<Ma.MQiwaiV)avt  a 


©  > 

NT  CD  < 

«/>  04  X 

X  ©  Ob  o»  O 

ob  ♦  O  o  * 

Or*  •  a 

o-  i  5  5 

•  2  □ 

r*  O  £  * 

w> 
a  r* 

O  v  *  •  * 

•  U  >  *  -J 

*  *  t  5  S 

U  H 

V)  >  VI  X  X  X  X 

g<U  3  «)  Irt  w 

zPo. 

«  vi  v>  vi 

t  «  H  »4  M 

!  sees 


AAAAAAAAAAAAi 
I  I  I  I  I  I  I  I  I  I  I  I 

i  i  i  »  i  i  i  i  i  i  i  i 

WiAr«»O^N','*rWiAf 


AAAAAAAAAAAAAAAAAAAAA 
»  I  I  •  I  I  «  I  I  I  I  t  I  I  I  »  I  1  I  I  I 
I  I  I  I  I  »  I  I  I  i  I  I  I  i  I  I  4  *  I  I  t 


AAAAAAAAAAAAA 
I  I  I  I  I  I  I  I  I  I  I  I  i 
I  I  I  I  I  I  I  I  I  I  I  I  I 


uuuuouuouvuuuuauouuulSuuuuuSuSuuLuuuuSuuuouulu 


64 


% 


o 


g 

3 


ill 


u 

in 

t 

5 


C 

a 

u 

*A 

3 

3 


i 


i  >  >  * 
~  o  o  0 


i  S  i 

3  3  3 

<A  iA  fA  i 

see 

V*  <A  v» 

A  A  A  A  A  A 

I  I  (  I  I  t 

I  I  I  I  I  I 


f 

a 


b  b 


2 

a 


ft 


«/> 

B  B 


8  £ 


A  A 

I  I 
I  I 


A  A 
I  I 
i  i 


.AAA 
I  t  I 
«  I  I 


AAAAAAAAAAAAAA. 
I  I  I  I  I  I  I  I  I  I  I  I  I  I 

I  I  4  I  I  I  I  I  I  I  I  I  I  t 


UOUL>OUUUUOOOOUUOU<JUOOOU<JOUOUUUUO< 


o  o  U  tj  u  o  t 


65 


TYPE :  120*4  CUN  TEST  CASE  CHAMBER  VOLUME  (Milt  0.00995  TRAVEL  <M)  I  4.75259 
GROOVE  DIAMETER  (Mil  0.11999  LAND  DIAMETER  (M)  t  0.1199*  CROOVE/  LAND  RATIO  (-)  »  1.000 
TNI  ST  (CALS/TURN):  *9.0  BORE  AREA  (M2):  0.01131  MEAT- LOSS  OPT  ION  t  1 
SHELL  THICKNESS  (Ml:  0.000114  SHELL  CP  (J/KC-KJ t  4(0.31(3  SHELL  DENSITY  (KC/M3) t  71(1.0920 
INITIAL  SHELL  TEMP  (K) :  2*3.  AIR  HO  (N2M«*2-K>:  11.3492 
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TTPE :  JA?  7P  CHAINS!  4803.1  1PF  HEIGHT  (XGI :  B.0000 
EROSIVE  COEFF  (-) :  0.000000  CHARGE  1GN  COOE:  0  CHARGE  IGN  AT  (S|t  O.OOOOOE.OO 
CRAIN  LENGTH  (Ml  s  0.016332  GRAIN  DIAMETER  (Ml!  0.009H4  PERT  DIAMETER  (Ml!  0.000308 
INNER  WEB  (Ml  :  0.001930  OUTER  NEB  (Ml:  0.001880 
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Ti  WITH  1  KJ/C  LLCCTHIC  1MVC?.S0«  DATE  TIME 
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B.3.  SAMPLE  CASE  3 


Here,  sample  case  2  has  been  modified  with  the  addition  of  propellant  bum  me  control.  The  bum 
rate  has  been  adjusted  to  produce  a  constant-breech-pressure  case.  The  gas  generation  rate  table  appeals 
in  the  SETC  deck  below,  after  the  electrical  power  input  deck. 

The  table  was  constructed  incrementally,  by  performing  a  series  of  runs,  and  adjusting  each  data  point, 
working  from  left  to  right,  to  smooth  out  the  pressure-versus-time  trace.  Note  that  this  could  not  have 
been  done  with  a  model  which  assumes  bum  rate  is  a  function  of  pressure  alone.  Over  the  long,  flat 
region,  the  gas  generation  rate  varies  from  800  to  4,000  kg/s,  even  though  the  mean  pressure,  on  which 
it  "depends,"  remains  nearly  constant. 


The  npu/  oeck: 


SC04S1 

anviirieo  ishvgj  icucwum  test  cut  « 

WITH  S  1  MS  rLAT  ELECTRIC  PULSE 

CONTAIN*  NO  A  TOTAL  Or  1KJ  ELECTA (C  PEA  OAAM  Or  PROPELLANT 
PROPELLANT  BURN  RAT"  TAILORED  TOR  fLAT  PRESSURE  PULSE 


siNro 

PUN  •  'T«,  .  KJ/O  ELEC-Ric.  DESIGNER  BURN  RATE'  OELT  -  JE-t  DELP  -  5E-3 


-NAD 

m 

2 

POPT  -  1.2.  1.0.  2 

SOFT  - 

0 

CfS 

m 

0.001 

•  HEAT 

T5HL 

9 

0.000114] 

CSHi.  •  460. 3143184 

A4HL  • 

7*41.0414 

TWAL 

9 

743 

HO  -  11 .34*318 

HL 

1 

SGUN 

NAME 

9 

'1J04O4  CUN 

TEST  CASK'  CHAM  -  0 

<• 

•> 

♦ 

o 

o 

MS  GAVE 

LAND 

9 

0.  .1 4444s 

G/L  -  1. 

TAAV  - 

4 .752544 

TNST 

9 

44 

SPAOJ 

NAME 

m 

'AfrSDS' 

PRIST  •  7 .0917  2 

SCOPS 

•POIS'  VALV.t  UStD  NITH  PARAMEDIC  PRINT  OPTION  POPTUl-J 
SPDI4 

CHON-' i  1AX  orCR-'OUT'  OIV-6494 .757 
IPO’S 

JhOW-'CH trr  OECK-'PROP'  NTH-J  0IV-0.4J354JJ7 

•  POIS 

SHOW-'j’AM'  OECK-'PROP'  NTH-J  DIV-O.CJ54 

•  POIS 

IHON-'PD'  DECK-' PROP'  NTH-J  DIV-O.OIJ4 
iPOIS 

SHOW-'  WES'  DECK-'  PROP'  NTH-J  OH/-..OJS4 

IPOIS 

3H0N-'  VMUZ'  TCK-'OUT'  0IV-0.304I 
SPOIS 

LH0  I?)'  OECK-'OUT' 

•Puts 

5H0N-'  LOUT  DECK- '  OUT' 

IRESI 

NP7S  •  4  AIR  •  I 

TIVRV  -  0,  O.OJOJJ,  0.074J,  4.744# 

PRES  •  O.SS44757,  17.J34M25,  0.4S447S7,  0.4S447S7 
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STDIS 

stdis 

STD  IS 

stdis 

STOIS 
STDIS 
STD  IS 
STDIS 
STDIS 
STDIS 
STDIS 


SHOW- 'TIME' 
SHOW-'TRAV' 
SHOW-' VEL' 
SHOW-'  ACCL' 
SHOW-' SUCH* 
SHOW-' MEAN' 
SHOW-'  BASE' 
SHOW-'EPWR’ 
SHOW-'  CEWE' 
SHOW-' III)' 
SHOW-' XBRN' 


OIV-.02S4  AEMK-' INCHES' 
MULT-3. 2801333  AEMX-'  FT/S' 
A1XX-' GRAVITIES' 
OIV-SIS4.7ST  RDK-'PSI' 
DIV-4894  .757  AIMX-'PSI' 
DIV-SSS4 .757  RIMR-'PSt' 
ROOC-' WATTS  ELECTRIC' 
AD4X-'  JOULES  ELECTRIC' 

RCMX-'  BURN  RATE' 


SACCO 

NAME  •  -  NONE' 


RECO  -  0  RCWT  -  0 


SPRIM 


NAME  - 

'RENITE' 

CHWT  -  0.001573944 

CAMA  - 

1.25 

rOAC  -  435174.7375 

COV  - 

0.001083819 

TEMP  -  2000 

NAME  • 

'  JA2  7P' 

CHWT  •  8. '3 

GRAN  • 

'  7pr* 

AHO  - 

1514. <11848 

CAMA  •  1.2257 

roue  • 

1142277.512 

COV  - 

NTSL-4 

0.000972778 

TEMP  •  3<00 

EROS  • 

0.0000000 

PR4L- 

13.789514,  27 

.579021,  48.94757, 

172.348925 

BA4L-  l 

0.02447,  0. 

038408,  0.074422. 

0.144424 

LEN  • 

0.0141322 

0IAM  •  0.010448 

P0 

•  0.000508 

WI  • 

0.0019104 

WO  •  0.0018794 

scorn 

OCriNC  A  TRAPEZOIDAL  ELECTRIC  POWER  PULSE 
100  HCOAWATTS  TOR  1  MILLISECOND. 

WITH  100  MICROSECOND  LINEAR  RAMPS  UP  AND  DOWN. 

SCTC 

NPWR  •  4 

PWR  -0.0.  »E9,  AES.  0.0 

TPWR-0. 000. 0.0001, 0.001. 0.0011 

NAAN-I 

TIRW-O. 000, 0.0001,0.0005,0. 001, 0.001S,  0.002,0.0025,  0.001, 
•AM-  10.0,  40.0,  200.0,500.0,1750.0,2500.0,1000.0,4500.0. 

SEND 


Produced  the  following  output: 


0.004 
5000. 0 
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SIMPLIFIED)  IMVC-2  BENCHMARK  TEST  CASE  ( 
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HJ/v.  tU-CTKIC.  DCS1GHCK  IUI«  GATE  lSHVCi.SO*  DATE  TIKE 
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KJ/C  ELECT*  1C.  DESIGNER  *U*N  RATE  IMVG2.6<M  DATE  TINE 
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B.4.  SK  E  CASE  4 


The  physics  of  this  case  is  identical  with  that  of  sample  case  3.  Changes  have  been  made  in  the 
STDIS  decks  to  display  the  trajectory  data  in  English  units.  The  purpose  is  to  illustrate  the  flexibility  of 
the  output  The  user  is  free  to  choose  the  units  of  his  choice. 

This  facility  is  provided  through  three  parameters:  MULT.  DIV,  and  REMK.  The  MULT  and  DIV 
parameters  are  constants  by  which  the  metric  variable  is  multiplied  or  divided  before  being  printed  in  the 
trajectory.  The  REMK  is  a  20-character  field  which  is  printed  a  the  top  of  each  page  next  to  the  column 
number  of  the  variable  being  displayed  in  the  trajectory.  This,  plus  the  12-character  name  of  the  variable, 
provides  enough  space  to  document  on  each  page  the  actual  units  in  which  the  variable  is  being  printed. 
It  is  of  course,  the  responsibility  of  the  user  to  insure  that  the  factor  is  consistent  with  the  claimed 
conversion! 


The  input  deck: 


scow 

SIMPLIFIED  IBHVC2  BENCHMARK  TEST  CARE  6 
WITH  A  1  M3  FLAT  ELECTRIC  PULSE 

CONTAINING  A  TOTAL  OT  1RJ  ELECTRIC  PER  CRAM  or  PROPELLANT 
PROPELLANT  BURN  RATE  TAILORED  TOR  PLAT  PRESSURE  PULSE 
TRAJECTORY  DATA  CONVERTED  TO  ENGLISH  UNITS 


S  INTO 


RUN 

• 

GRAD 

- 

CPS 

• 

SHCAT 

TSHL 

- 

TOAL 

- 

SCUN 

name 

• 

LAND 

» 

TWIT 

m 

SPROJ 

NAME 

m 

SRXSI 

NPTS 

- 

TRAV 

• 

PRES 

• 

STDIS 

SHON-'TIMC' 

STDIS 

SHOW-'TRAV' 

STOIS 

SHOW-' VEL' 
STDIS 

SHOW-'  ACCL' 
STDIS 

SHON-' 9RCH' 
STOIS 

SHOP-' MEAN' 
STOIS 

SHOW*' BASE' 
STOIS 

SHOW-'EPWR' 


'T6,  1  KJ/C  ELECTRIC,  OESICNER  BURN  RATE'  OCLT  -  SE-6  DELP  -  5E-5 
2  POPT  -  1,2, 1,0. 2  SOFT  •  0 

0.001 

0.0001143  CSHL  •  460.3163186  RSHL  •  7861.0916 

293  HO  •  11.348218  HL  •  1 

'  120404  GUN  TEST  CASE'  CHAM  •  0.009946948  GRVE  •  0.1199896 
3.1199896  C/L  -  1.  TRAV  •  4.752S94 

99 

■APrSDB'  PRNT  •  7.09872 

4  AIR  -  1 

0,  0.02032,  0.0762,  4.7498 

0,6894757,  17.2368925,  0.6894757,  0.6694757 

REMK*'  MILLISECONDS' 

DIV-0.02S4  REMX-* INCHES' 

MULT-3. 2808333  REMK-' FT/S' 

DIV-9.8067  REMK-' GRAVITIES' 

DIV-6894.7S7  REMK-' PRESSURE, PSI' 

OIV-6894 .757  REMK-' PRESSURE,  PSI' 

OIV-6894.7J7  REMK-' PRESSURE,  PSI' 

DIV-746.0  REMK-' HORSEPOWER,  ELECTRIC' 
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STD  IS 
STDIS 
CTDIS 
SACCO 
SPAIM 

STDIS 

STOIS 

STDIS 

STDIS 

STDIS 

JRCCO 


SHOW-'EENC'  MULT-1. 356  RXMK-' TOOT-POUNDS  ELECTRIC' 
SHOW-' Z ( 1 J  * 

SHOW—' XBRN*  DIV-0.0254  REMK- ' BURN  RATE  (IN/SEC)' 
NAME-' NOME'  PECO  •  0  RCWT  •  O 


NAME  •  'SEMITE' 

GAMA  -  1.25 

COV  •  0.001083819 


CHWT  -  0.0015739*6 
TORC  -  83517*. 7375 
TEMP  -  2000 


SHOW-'  BASE' 
SHOW-'  EPWR' 
SHOW-'  ECMt' 
SHOW-' 2(1)' 
SHOW-'  XBRN' 


0IV-8894 .757  RtMX-'PSI' 
RESOt-' WATTS  ELECTRIC' 
RIMR-' JOULES  ELECTRIC' 

REMR-'SURM  RATE' 


NAKZ 


'HOME' 


RICO  -  O 


RCWT  -  0 


SSPRIM 

NAME  •  '  BENITE'  CHWT  •  0.00157396* 

CAMA  -  1.25  fORC  -  63517*. 7375 

COV  •  0.001083819  TCMP  •  2000 

STROP 

NAME  -  '  JA2  7P'  CHWT  -  8.0  CRAM  -  '  7PP' 

RHO  •  1586.611868  CAMA  -  1.2257  COAC  •  1142277.932 

COV  •  0.000992778  TOW  •  3400  EROS  -  0.0000000 

NTBL-4 

PR4L-  13.789514.  27.579028,  48.94757,  172.368925 

0R4L-  0.02667,  0.038608,  0.074422,  0.144624 

LCN  -  0.01*3322  OIAM  -  0.010**8  PD  -  0.000508 

MI  •  0.0019304  WO  •  0.0018794 


SC0991 

DEFINE  A  TRAPEZOIDAL  ELECTRIC  POWER  PULSE 
100  MEGAWATTS  TOR  1  MILLISECOND, 

WITH  100  MICROSECOND  LINEAR  RAMPS  UP  AND  DOWN. 

SCTC 

NPWR  -  4 

PWR  -0.0,  8C9.  IE9.  0.0 

TPWR-0. 000, 0.0001, 0.001, 0.0011 

N8RM-9 

T8RN-0. 000.0. 0001, 0.0005,0. 001, 0.0015.  0.002,0.0025,  0.003, 
8RM-  10.0,  60.0,  200.0.500.0.1750.0.2500.0,3000.0,4500.0, 

SEND 


Produced  the  following  output: 


0.004 

5000.0 
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TYPE :  120!*1  CUN  TEST  CASE  CHAMBER  VOLUME  (Mil:  0.0CJ3S  TRAVEL  (HI)  :  4.7S2S* 
CROOVE  DIAMETER  (Ml  :  0.I19»«  LAND  DI NHETEA  (H)  :  0.11»»»  CROOVC/LAND  RATIO  (-)  :  1.000 
TWIST  (CALS/TURN):  »».0  BORE  ARCA  (M2):  0.01131  HEAT-LOSS  OPTION:  1 
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BURN  INC  RATES  (M/S) :  -  -  -  0.02(4? 

BURNING  RATES  (M/S) ;  -  -  -  0.03141 

BURNING  RATES  (M/S) t  -  -  -  0.01442 

BURNING  RATES  (M/S) :  -  -  -  0.14442 


fHA JCCTORY  VARIABLES:  /  1/  TRAJ  1  TIME 


V)  h 
uj  — 

X  tA  > 
U  'v  < 

Z  f>a  i 
-  U  O  t 


G 

a 

ui 

*A 

*■**-**-*£ 


>  Jzzumu-. 
HXflZsuyNB 


232255252:2 


NNSNNSVSVm 


O  «/*  O  OOOOOOOOOOOOOO'-t  O  OOOOOOO  OOOO®  oo 

O  ®  O  OOO^nO'^Oi/^OOOOOO'*  o  OOOOOOO  .  ■  . 

O  «*»  O  O  >n»  O  •  .......  o  "fl  O  *ni  c*>  0«/> 

. ilNOOOOO  Orfi  o  0  0  0  0*^10^1  o  NtAr*#  On 

o  r~  >/>  o  c*  a>  o  <n  «/*  o  ^  n*  m  .z'  ® 

— •  —*  — ''n  ^  ^  */>  s«  i“*  ®  .-«  --«  *-> 

•  l"»  nmNNNNNNMWNX*****  —  ^  *-  mm  <—  _  ^  ^ 

O  O  OOOOOOOOOOOOOOO  O  OOOOOOO  ooooe  oo 

I  I  I  I  I  I  I  i  I  I  t  (  I  I  I  I  I  I  1  t  I  I  I  I  I  I  I  I  I  I  11 

ui  u  uuuubiuuuuuyuuuu  ui  uuiyiJuuu  u  u  ui  u  ui  u  u 

O  ^  n  o^Nr*i/»r,‘Nr(inrt(#i-p«*iO  n  w  •*  ft  w  */»  **i  nc 

am  >#  — •  «nttr>NNO  *ONNn  «-•  r~ 

O  ®  N»r»*n^r*  a>  */> 

o  xr  r»  on^cnyH  p- ♦-*«*,*  cs  • 

0«n  •-*  «  r«NNN(*m»  «<n««r>  r- ® 

OO  O  OOOOOOOOOOOOOOO  O  OOOOOOO  ooood  oo 

%•  n*  »•  c»  r»  r»  r- r»  r»  r»  r*  r»  f*  *»•  f*  f»  r*  r  f*  f*  f*  f*  ^  r-  r*  r-.  r»  r-  r*  r* 

O  OOOOOOOOOOOOOOO  O  OOOOOOO  ooooo  oo 

♦  ♦  ♦  ♦♦♦♦♦♦♦  »  •  «  «  •  ♦  * 

UI  UlUIUlUIUlUlUUIUlUlUIUIUlUlUl  UI  UhlUUUUU  UI  UI  UJ  UI  UI  UI  UI 

o  f*  O  0  0  0  0  0  0  0  0  0  0  0  0  0  0--  O  OOOOOOO  OOOOO  OO 

a  *  o  ooooooooooogoowi  a  ooooooo  ooooo  oo 

2**  2  goooooooogogoof-  o  ooooooo  ooooo  oo 

or*  O  OON«O*0N>«Ov«N<««  o  «aiM4»oo  ooooo  oo 

O  y  «  «  s0  <Ar*>  r  r*  a  O  OSSSO  <x>  00 

•  Off  •  *•••••»•«•••»•«  ♦  I  I  I  ■  f  ■  • 

o  O  OOOOOOOOOOOOOOO  o  ooooooo  OOOOO  Oo 

O  O  OOOOOOOOOOOOOOO  o  ooooo 

♦  *  ♦♦*♦♦♦♦♦♦♦♦♦*♦♦  *  »  »  «  «  • 

ui  u  u  uuuuuy  uguuuuu  u  ui  guuuu 

Ow  O  OOOOOOOOOOOOOOO  o  ooooooo  ooooo  oo 

oa  O  OOOOOOOOOOOOOOO  O  ooooooo  ooooo  oo 

o  r>  o  OOOOOOOOOOOOOOO  O  OOOOOOO  ooooo  oo 

O  —  O  OOOOOOOOOOOOOOO  o  ooooooo  ooooo  oo 

O  *-*  v  (SO0OogDll)(SID«aiOMSB<S  CD  SOQOvOO  OOOOO  oo 

OO  O  OOOOOOOOOOOOOOO  o  OOOOOOO  OOOOO  oo 

m  .  .......  . 

f*  O  n  •  •  •  *—  o»  ♦  n  « i/i  o  m  >o  ^  <d  o*  a> 

ifl  *  <**o«cNf-O*^*#f*-*0«*>r,,*DOt/»  r-  /li^  v  h0i  n  o  o  r<  as  ^  n  v 

•  O  r*  m»*inior*®o»o^n<iNNNtfl  O  •"I'-OOO  a  —■ 

m  o  — •  N«<cr  vac  od  <d  cd  eo  <s>  a  c 

nn  **  «/l  —  «n  nn  n  n 


— •  n*  •* 


o  — • 
O  o 
«  r* 
<-x  *n 


* 

«n 
«n  ^ 


O  IM 
o  *# 
*v  o 


<noo<n^Or 


>C*H«n*D<A<A  *-* 


>  «st  o  o 

4  m  r>l 

n»«0 
r  ^  ^ 


Nmo 
NO® 
•  O  *D 

«-  O*  — 
tn«n  in 


mnni/iin 


o  o 

r*>  w» 
cn  *a 


w  _ 

^  m 


'*>0v0/iD®r-v«40«  «-*  rx®nn 

w  <-»  t*  m  «  Mi  -*  m  O  *•  «  r-  O  kin^n 

r  r*-  O*  cm  d>  *r»  r-  r*  r-  ex  «o  nflin-fl 

r»  o  omop* 

-iM«N(XNnrw*t*  v  *n  */>  *n 


®  i/i  r* 
—  «**./> 


i.i  >n  in  n  in 


o  o 
o  o 
o  o 
o  o 
o  a 

o  a 


o  o 
o  o 
o  o 
oo 
o  o  t 

o  o  I 


o  o 
o  o 
o  o 
o  o 
o  o 


,2 


fOn*  •  •  i  •  . 

■  •  •  ■0'"*«3fcW»*D0r*000® 

^  i*  r~  vi/i® 

Oni/noifi®ON'^i^r<or>  -4  *n 
Hi«  ®ON*  ^r«  OO 


▼O’*  O  r* 
*  K) 

«ON«®  f  > 

y  r-  O  ^  /M  w  /*, 

m  ♦  v  «r  • 


«00«N 
•n  —  rsi  00  CD 
n»  vO  «o  w 

4NNNM 


vD  vAl 

*n 
♦  «/l 
nx  «N 


r>  (t  4  in  y  «r*  iMo^NnovtfNv 
iifti**Nn0^n  *  •  •  •  ♦  •  •  • 
m  rt  A  •  .  •  •  tftOBNNinNM 

N  I  •  ^Ohr>«**n®Niioifl® 


v  O  ®  n  in 


OL  nx 
D  f* 

<0  rx 


JOOOOOO 
E  I  »  I  t  I  I 
J  (J  UI  (J  UI  UI  UI 

»NNNin®inpx®ino^-"Nnin 

}<s^v'*i«#BiDin«o<«SMCoa 

14^0  aiNrN0iNr0i*«OO 


G 

a 

UI 

O  + 
X  ^ 


'3'"" 


<N**Tr*-4«a 

2  ■  r . 


^  ^  o  r*  ^ 

*#  r*  ®  2 


*#*n  r>  *d  o  < 
Ono4f  * 
N«Xi 
<  ' 


'  •  n  f«n  n  o  m 

iininc#!-  v/i 

*  •  O  * 


■3  — 

C 

ui 

h 


i  ®  D-  rv  O  O  O  ® 
w  *4»  ®  «n 

-*Nr  ®  7r*o 
1 in  ix  k  s 


cooooodoooodoo 


«  ^»«NNNWf»  ^  ^  in  XI  X)  B 


o  o 

'v  on 

—  o  r~ 


»-  i 

a  uj 

<  o 

h*  o 

CO  O 

'  I  o 

*  f-  i r> 

-  O  * 

>  X  o 

</> 


aoooooooooooooo©' 

OOOOOOOOOOOOOOO 

-loaooooooooooooc* 

uononononononon/i 

l**«Nfwnn#x^i/ii(xnM* 

a: . . . 

-eooooooooooooooo  « 


X  O  c 
a.  O  c 
O 

•J  o  . 


<  D  <c  a>  ®  ® 


•  ->  o  o  o  o  a. 

ooooo  a. 
’OOnon 
)noo^< 


looonion 
>  O  O  O  0.00 
(nonn  o 


- - ^ - ‘ 


89 


90 


3 . 2000  59.350 


TKAJIOTOHT  VARIABLES :  /  1/  TRAJ  1  TIME 


OOOCOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOO 

. ooooooooooooooooooooooooooo 

^0^0>A00  OOOOOOOOOOOOOOOOOOOOOOOOOOO 

r-0<V^r>00  OOOOOOOOOOOOOOOOOOOOOOOOOOO 

v®  r-  r~  r>  r*  <s  O  OOOOOOOOOOOOOOOOOOOOOOOOOOO 

O  OOOOOOOOOOOOOOOOOOOOOOOOOOO 


V  r*  r>*  r* 


m  ~  *  pm  O  •  O 

mmr*  r-  r+  o 


ooooooooooooooooooooooooooo 

ooooooooooooooooooooooooooo 

- oooooo 

o  o  o  o  o  o 


OOOOOOOOOQOOOOQOOQOOO 

OOOOOOOOOOOOOOOOOOOOO 


OOO  oooooo 


OQO  OOOOOOO  QOOOOO< 


>0000000000000000000 


uuu 

OOO 

OOO 

O  o  o 

OOO 

ODtt« 


UUhlUUUU 

OOOOOOO 

OOOOOOO 

OOOOOOO 

OOOOOOO 

8SAC90O 


UUUuJUUUb!  UUhlUUUUUUbiUUUUUUUUU 
OOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOO 


OOO  OOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOO 


OOO 

OOO 

OOO 

OOO 

OOO 


OOOOOOO 

OOOOOOO 

OOOOOOO 

OOOOOOO 

OOOOOOO 


OOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOO 


OOO  OOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOO 


H 

< 

o 


®Q« 


®«r  r* 


ui  — 

X  IA  > 

U  M.  <  ~ 

ZHCini 

-Luai 


u  — 

-•  ac 

O  Ui 
UI  -i 

-J  u 

UJ 

«/> 


5 

o 


>n«is 
r-  00 


NNNrvMMfR 
«/>«/>  i/>  i/> 


p*  (M  r*  N  ®  r* 

«*»»»»» 

m  ift  </)  «n  «n  irt 


fM*®0^s®r*r-‘p^r*««/>®r 

r*~.SkP'’tPMO«'#tf><"tPM-«*fl 


>  JZZUIItU^Z 

U(jzuUa,U''S) 

hxaiauuNx 

SS222SS222 

Hf-fr-Hf-HHHHH 


SVSN' 


.  NVSSNSN 


}  n  o  -h  /i  rs 

»  O  O  ®  r*  ▼  ®  fv  OJ  py 

>  O  O  O  ^  a>  r*  a>  ao 

p' pi  (v  ry  <y  M  — 


Mfl®  #®®r  MVtf'P'rflP'  0»*»'-®N0/'®®><PS»/l''r'  ®  V 


O®  C  P*  ®  ®  vOAiAV^P'O^rbO*® 

p*nor» 

«®00«NN 


*®N®OV®NiO(>N/'^P»in 

®»OO"',«0((VNnP'p'*» 


H 
UI 

r*  fl®N®  <VO® 


■^0D^/i<-QO»rv*9>'nO'*>P>'pP'r 
>®Mf''\«/'UfflnO®P’^<Pf'  •  •  • 


l^aop'/i®  ^  1 _ 

•®  v®  r~  r—  p-  p-*  r~  3®aaa^®®9  • 
J  ® 


O  O  pm 
O  O  (T* 
V*  O  PM 


a.  O  O  O  O  O  O  r 


tOOOOOOOOOO 


OOOOOOOOOOOOOOOO 
-OOOOOOO  -JO  000000000  -oooooooooooooaoo 
■ziO/'O'/’OPPJ'^O^O'P’O/’O^O^O'P'O/'O^O'/iOi/'O'P'O'/'O'^ 
J  n  y  *  i/i  1/1  ifl  j)  iJ  so  P“-  p-  asilM>00’4MNN^P)vVi/<</>4>AP  r~  ®  O  ®  ® 


91 


TRAJECTORY  VARIABLES 
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TIME  (MSI:  3.329  5.383 
TRAVEL  (Ml:  1.4552  4. 1524 
VELOCITY  (M/SI  1170.84  1732.84 
ACCELERATION  (C| :  43100.  14899. 
BREECH  PRESS  (MPA):  410.4479  148.2100 
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No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  Commander 

U  S.  Army  Materiel  Command 
ATTN:  AMCAM 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

1  Director 

U  S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-AD, 

Tech  Publishing 
2800  Powder  Mill  Rd. 

Adelphi.  MD  20783-1145 

1  Director 

U  S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-O-AD, 

Records  Management 
2800  Powder  Mill  Rd. 

Adelphi.  MD  20783-1145 

2  Commander 

U.S.  Army  Armament  Research. 

Development,  and  Engineering  Center 
ATTN:  SMCAR-IMI-1 
Picatinny  Arsenal.  NJ  07806-5000 

2  Commander 

U.S.  Army  Armament  Research. 

Development,  ard  Engineering  Center 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
U.S.  Army  Armament  Resed/ch, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCB-TL 
WatervUet,  NY  121894050 

1  Director 

U.S.  Army  Advanced  Systems  Research 
and  Analysis  Office  (ATCOM) 

ATTN:  AMSAT-R-NR.  MJS  219-1 
Ames  Research  Center 
Moffett  Field.  CA  94035-1000 


1  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMJ-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

1  Commander 

U.S.  Army  Tank- Automotive  Command 
ATTN:  AMSTA-JSK  (Armor  Eng.  Br.) 
Warren,  MI  48397-5000 

1  Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSR 

White  Sands  Missile  Range,  NM  88002-5502 

(Om*.  only))  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 

Fort  Benning,  GA  31905-5660 


(Und**.  only)  i  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-WCB-O 
Fort  Benning,  GA  31905-5000 


1  WL/MNOI 

EglinAFB.  FL  32542-5000 

Atenton.  Provinfc£rcupd 

2  Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

1  Cdr,  USATECOM 
ATTN:  AMSTE-TC 

1  Dir,  USAERDEC 
ATTN:  SCBRD-RT 

1  Cdr,  USACBDCOM 
ATTN:  AMSCB-CB 

1  Dir,  USARL 

ATTN:  AMSRL-SL-I 

5  Dir,  USARL 

ATTN:  AMSRL-OP-C1-B  (Tech  Lib) 
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No,  of 

Copies  OrcaniEatiqn 

1  OSD/SDIO/1ST 

ATTN:  Dr.  Len  Caveny 
Pentagon 

Washington,  DC  20301-7100 

1  Director 

U.S.  Army  BMD 
Advanced  Technology  Center 
P.O.  Box  1500 
Huntsville.  AL  35807 

1  Chairman 

DOD  Explosives  Safety  Board 
Room  856-C 
Hoffman  Bldg.  1 
2461  Eisenhower  Ave. 

Alexandria,  V A  22331-0600 

1  Deportment  of  the  Army 

Office  of  the  Product  Manager 
15Smm  Howitzer,  M109A6.  Paladin 
ATTN:  SFAE-AR-HDMP,  Mr.  R.  De  Kleine 
Picatinny  Arsenal,  NJ  07806-5000 

1  Commander 

Production  Base  Modernization  Agency 
U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  AMSMC-PBM-E,  L.  Laibson 
Picatinny  Arsenal,  Ni  07806-5000 

1  PEO-Armaments 

Project  Manager 
Tank  Main  Armament  Systems 
ATTN:  AMCPM-TMA,  K.  Russell 
Picatinny  Arsenal,  NJ  07806-5000 

1  PEO- Arman,  nts 

Project  Manager 
Tank  Main  Armament  Systems 
ATTN:  AMCPM-TMA- 105 
Picatinny  Arsenal,  NJ  07806-5000 

1  PEO-Armamcnis 

Project  Manager 
Tank  Main  Armament  Systems 
ATTN:  AMCPM-TMA- 120 
Picatinny  Arsenal,  NJ  07806-5000 


No.  of 

Copies  Organization 

6  Director 

Benet  Laboratories 
U.S.  Army  Watervliet  Arsenal 
ATTN:  SARWV-RD 
L.  Johnson 
G.  Carafano 
R.  Thierry 
R.  Hasoenbetn 
P.  Votis 
P.  Alto 

Watervliet,  NY  12189 

3  Commander 

U  S.  Army  AMCCOM 
ATTN:  AMSMC-1RC,  G.  Cowan 
SMCAR-ESM(R), 

W.  Fortune 
R.  Zastrow 

Rock  Island,  IL  61299-7300 

1  Commander,  USACECOM 
R&D  Technical  Library 
ATTN:  ASQNC-ELC-IS-L-R,  Myer  Center 
Fort  Monmouth,  NJ  07703-5301 

1  Commandant 

U.S.  Army  Aviation  School 
ATTN:  Aviation  Agency 
Fort  Rucker,  AL  36360 

1  Director 

HQ,  TRAC  RPD 

ATTN:  ATCD-MA,  MAJ  Williams 
Fort  Monroe,  V A  23651-5143 

1  Headquarters 

U.S.  Army  Materiel  Command 
ATTN:  AMCICP-AD,  Michael  F.  Fisette 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

1  Commander 

U.S.  Army  Armament  Research 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCS 
Picatinny  Arsenal,  NJ  07806-5000 
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No.  of 

Copies  Organization 

1  Commander 

U.S.  Army  Armament  Research 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCH-T,  L.  Rosendorf 
Picatinny  Arsenal,  NJ  07806-5000 

1  Commander 

U.S.  Army  Armament  Research 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCH-V,  E.  Fennell 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

U.S.  Army  Research  Office 
ATTN:  Technical  Library 
D.  Mann 
P.O.  Box  1221 

Research  Triangle  Park.  NC  27709-2211 

1  Commander 

U.S.  Army  TRAC  -  Fort  Lee 
Defense  Logistics  Studies 
Fort  Lee.  VA  23801-6140 

1  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AE.  J.  Picard 
Picatinny  Arcenal,  NJ  07806-5000 

5  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
aTTN:  SMCAR-AEE-B, 

A.  Beardel! 

D.  Downs 

S.  Einstein 
A.  Bracuti 
D.  Chiu 

Picatinny  Arsenal,  NJ  07806-5000 

I  Commander 

U.S.  Army  Armam:ni  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AEE,  /.  Lannon 
Picatinny  Arsenal,  NJ  07806-5000 


No.  of 

Cacits  fltnaaiaiifla 

1  Commander 

U.S.  Army  Armament  Research. 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AES,  S.  Kaplowitz 
Picatinny  Arsenal,  NJ  07806-5000 

1  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-FSA-T,  M.  Sa'sbury 
Picatinny  Arsenal.  NJ  07806-5000 

6  Commander 

U.S.  Army  Armament  Research. 

Development,  and  Engineering  Center 
ATTN:  SMCAR-FSE, 

T.  Cora 

B.  Kiiutelsky 
K.C  .  Pan 
W.  Davis 

C.  Durham 
A.  GraT 

Picatinny  Arsenal.  NJ  07806-5000 

6  Commander 

U.S.  Army  Armament  Research. 

Development,  and  Engineering  Center 
ATTN:  SMCAR-EG. 

G.  Ferdinand 

H.  Naber-Libby 
R.  Lund  berg 

J.  NiJes 
R.  Moreira 
W.  Morelli 

Picatinny  Arsenal.  NJ  07806-5000 

1  Commander 

U.S.  Army  Belvoir  R&D  Center 
ATTN:  STRBE-WC 

Technical  Library  (Vrult) 

Bldg.  315 

Fort  Belvoir.  VA  22060-5606 

1  President 

U.S,  Army  Artillery  Board 
Fort  Sill,  OK  73503 
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No.  of 

£2Qi£3  Organization 
1  Commandant 

U.S.  Army  Command  and  General  Staff 
College 

Fort  Leavenworth,  KS  66027-5200 

1  Commandant 

U.S.  Army  Special  Warfare  School 
ATTN:  Rev  and  Tog  Lit  Div 
Fort  Bragg.  NC  28307 

1  Commander 

Radford  Army  Ammunition  Plant 
ATTN:  SMCRA-QA/HI  Library 
Radford.  VA  24141 

1  Commandant 

U.S.  Army  Field  Artillery  School 
ATTN:  STSF-TSM-CN 
Fort  Sill  OK  73503-5600 

4  Deputy  Commander 

Strategic  Defense  Command 
ATTN:  SFAE-SD-HVL, 

S.  Smith 
LTC  Kee 
D.  Lianos 

T.  Aden 
P.O.  Box  1500 
Huntsville,  AL  35887-8801 

2  Commander 

U.S.  Army  Foreign  Science  and  Technology 
Center 

ATTN:  AMXST-MC-3. 

S.  LeBeau 
C.  Beiter 
220  Seventh  St,  NE 
Charlottesville,  VA  22901 

1  Commandant 

U.S.  Army  Field  Artillery  Center  and  School 
ATTN:  ATSF-CO-MW,  B.  Willis 
Fort  Sill.  OK  73503 

I  Naval  Sea  S).:t».m  Command 
Depafjneri  of  the  Navy 
ATIN:  CSEA,  CDR  Dampier 
06KR12 

Washington,  DC  20361-5101 


No.  of 
Copies 

1 


2 


1 


1 


2 


2 


3 


1 


Organization 

Office  of  Naval  Research 
ATTN:  Code  473,  R.  S.  Miller 
800  N.  Quincy  St 
Arlington,  VA  22217 

Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA62R 
SEA  64 

Washington,  DC  20361-5101 
Commander 

Naval  Air  Systems  Command 
ATTN:  AlR-954-Technical  Library 
Washington.  DC  20360 

Naval  Research  Laboratory 
Technical  Library 
Washington,  DC  20375 

Commander 

Naval  Surface  Warfare  Center 
ATTN:  J.  P.  Consaga 
C.  Gotzmer 

Silver  Spring,  MD  20902-5000 
Commando 

Naval  Surface  Warfare  Center 
ATTN:  K.  Kim/CodeR-13 

R.  Bemecker/Code  R-13 
Silver  Spring.  MD  20902-5000 

Commander 

Naval  Surface  Warfare  Center 
Indian  Head  Division 
ATTN:  610,  C.  Smith 
61IQJ,  K.  Rice 
61  IOC,  S.  Peters 
Indian  Head,  MD  20640-5035 

Commander 

Naval  Underwater  Systems  Center 
Energy  Conversion  Dept 
ATTN:  Technical  Library 
Newport  RI  02840 
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No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


6  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  G33. 

T.  Doran 
J.  Copley 

Code  G30,  Guns  and  Munitions 
Division 

Code  G301,  D.  Wilson 
Code  G32.  Gun  Systems  Division 
Code  E23,  Technical  Library 
Dahlgren,  VA  22448-5000 

1  Commander 

Naval  Weapons  Center 
ATTN:  Code  388,  C.  F.  Price 
Info  Science  Div 
China  Lake,  CA  93555-6001 

1  OLAC  PL/TSTL 

ATTN:  D.  Shipleu 
Edwards  AFB,  CA  93523-5000 

10  Central  Intelligence  Agency 

Office  of  Central  Reference  Dissemination 
Branch 

Room  GE-47  HQS 
Washington.  DC  20502 

1  Central  Intelligence  Agency 
ATTN:  Joseph  E.  Backofen 
NHB,  Room  5N01 
Washington,  DC  20505 

5  Director 

Sandia  National  Lzooratories 
ATTN:  T.  Hitchcock 

R.  Woodfin 
D.  Benson 

S.  Kempka 
R.  Beasley 

Advanced  Projects  Div  14 
Organization  9123 
Albuquerque,  NM  87185 

2  Director 

Los  Alamos  National  Laboratory 
ATTN:  B.  Kaswhia 
H.  Davis 

Los  Alamos.  NM  87545 


1  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  M.  S.  L-355,  A.  Buckingham 
P.O.  Box  808 
Livermore.  CA  94550 

2  Director 

Sandia  National  Laboratories 
Combustion  Research  Facility 
ATTN:  R.  Armstrong 
S.  Voaen 
Division  8357 

Livermore,  CA  94551-0469 

i  University  of  Illinois 

Dept  of  MechyindusL  Engr. 

ATTN:  Prof.  Herman  Krtcr,  144  MEB 
1206  N.  Green  St 
UrUana.IL  61801 

1  Johns  Hopkins  University /CP  LA 

ATTN:  T.  Christian 
10630  Little  Patuxent  Parkway,  Suite  202 
Columbia,  MD  21044-3200 

1  Pennsylvania  State  University 

Dept,  of  Mechanical  Engr. 

ATTN:  Dr.  K.  Kuo 

312  Mechanical  Engineering  Bldg. 

University  Park,  PA  16802 

1  North  Carolina  State  University 

ATTN:  John  G.  Gilligan 
Box  7909 

1 1 10  Burlington  Engineering  Labs 
Raleigh,  NC  27695-7909 

1  State  University  of  New  York  at  Buffalo 

Department  of  Electrical  Engineering 
ATTN:  Dr.  W.  J,  Sargeant 

Bonner  Hal!  •  Room  312 
Buffalo,  NY  14260 

2  institute  for  Advanced  Technology 
ATTN:  Dr.  H.  Fair 

Dr.  T.  Kichne 
4030-2  West  Baker  Lane 
Austin,  TX  78759-5329 


99 


No.  of 

Copies  Organization 

1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Technical  Library 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 

1  SPARTA 

ATTN:  Dr.  Michael  Holland 
9455  Towne  Center  Dr. 

San  Diego.  CA  92121-1964 

5  FMC  Corporation 

ATTN:  Mr.  G.  Johnson 
Mr.  M.  Seale 
Dr.  A.  Giovaneui 
Mr.  J.  Dyvik 
Dr.  D.  Cook 
4800  East  River  Rd. 
Minneapolis,  MN  55421-1498 

10  GT  Devices 

ATTN:  Dr.  S.  Goldstein 
Dr.  R.  J.  Greig 
Dr.  N.  Winsor  (8  cp) 
5705A  General  Washington  Dr. 
Alexandria.  VA  22312 

3  General  Dynamics  Land  Systems 
ATTN:  Dr.  B.  VanDeusen 
Mr.  F.  Lunsford 
Dr.  M.  Weidner 
P.O.  Box  2074 
Warren.  MI  48090-2074 

2  Alliant  Techsymcms,  Inc. 

ATTN:  R.  E.  Tompkins 

J.  Kennedy 
MN38-330 

10400  Yellow  Circle  Dr. 
Minnetonka.  MN  55343 

2  Olin  Ordnance 

ATTN:  V.  McDonald,  Library 
Hugh  McElroy 
P.O.  Box  222 
c».  Marks.  FL  32355 


No.  of 

Copies  Organization 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.  S.  Gough 
1048  South  Sl 

Portsmouth.  NH  03801-5423 

1  Physics  International  Library 
ATTN:  H.  Wayne  Wampler 
P.O.  Box  5010 

San  Leandro.  CA  94577-0599 

2  Rockwell  International 
Rocketdyne  Division 
ATTN:  BA08, 

J.  E.  Flanagan 
J.  Gray 

6633  Canoga  Ave. 

Canoga  Park.  CA  91304 

2  Princeton  Combustion  Rsch  Lab 

ATTN:  M.  Summerfield 
N.  Messina 

Princeton  Corporate  Plaza 
1 1  Dccrpark  Drive 
Bldg.  IV.  Suite  119 
Monmouth  Junction.  NJ  08852 

2  Science  Applications.  Inc. 

ATTN:  J.  Batteh 

L.  Thornhill 
1519  Johnson  Feny  Rd. 

Suite  300 

Marietta.  GA  30062-6438 

I  Eli  Freedman  &  Associates 

ATTN:  E.  Freedman 
2411  Diana  RdL 
Baltimore  MD  21209 

I  Veritay  Technology,  Inc. 

4845  Millenport  Hwy. 

P.O.  Box  305 

East  Amherst,  NY  14051-0305 

1  Baitclle 

TWSTIAC 
505  King  Ave. 

Columbus,  OH  43201-2693 
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2  California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
ATTN:  L.  D.  Strand.  MS  123-224 
D.  Elliot 

4300  Oak  Grove  Dr. 

Pasadena,  CA  91109 

1  Oenerai  Electric  Co. 

Defense  Systems  Division 
ATTN:  Dr.  J.  Mandzy 
Mail  Drop  43-220 

100  Plastics  A  ve. 

Pittsfield.  MA  01201 

2  SAIC 

ATTN:  Mr.  N.  Sinha 
Dr.  S.  Dash 
301  Office  Center  Dr. 

Fort  Washington.  PA  19034*3211 

1  Science  Application  International  Corporation 
ATTN:  Dr.  O.  Chryssamellis 
8400  Nomundele  Blvd. 

Suite  939 

Minneapolis.  MN  33437 

Atafecn  -Proving  Ground 

4  Commander,  USACSTA 
ATTN:  S.  Walton 
G.  Rice 
D.  Lacey 
C.  Herud 
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2  RARDE 
GS2  Division 
Bldg.  R3I 

ATTN:  Dr.  C.  Woodley 
Dr.  G.  Cook 
Port  Halstead 

Sevenoaks,  Kent  TNI 4  7BP 
England 

2  Director 

Fraunhofer  -  [nation  far  Kurzzeiidynamik 
Ernst  -  Mach  Ins  ti  tut 
ATTN:  Dr.  G.  Sc  breeder 

Prof.  G.  Klingcrberg 
Hauptstrafie  18 
D-7858  Weil  am  Rhein 
Germany 

3  Soreq  Nuclear  Research  Center 
Israel  Atomic  Energy  Commission 
ATTN:  Dr.  D.  Saphier 

Dr.  R.  Alimi 
Ms.  D.  Plotnik 
Yavne  70660.  Israel 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-348 _ Date  of  Report  January  1994 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS _ 

Street  or  P.O.  Box  No. 


City,  Stale,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 
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BUSINESS  REPLY  MAIL 

OUT  CUSS  KMt  ta  tan.  vs.  to 


Po*uq«  «ill  b«  and  Oy  addrassM 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Ub) 
Aberdeen  Proving  Ground,  MD  21005-5066 


NO  PCSTACt 
MCISSAAV 
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